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Notations and conventions 



a lattice spacing 

A }1 = A a p t a gauge field 

(A 2 ) A 2 -condensate 

olf («g) infrared exponent of the ghost (gluon) propagator 

jS = ^ bare lattice coupling 

So 

P(g), Pi the renormalisation group beta function and 

its first coefficients 

lis) ' 7> 7i anomalous dimension ; anomalous dimension 

in a generic MOM scheme and its coefficients 

e p a vector in direction ]i of norm a 
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General introduction 



Of course, you can put a theory on the lattice. But then - it is a mess !" 

Giorgio Parisi, les Houches 



This PhD dissertation is devoted to a non-perturbative study of QCD correlators 
in Landau gauge. The main tool that we use is lattice QCD. It allows a numerical 
evaluation of the functional integrals defining vacuum expectation averages of the 
theory i.e. Green functions. The advantage of this method is that it gives access to 
the non-perturbative domain and exactly preserves the gauge-invariance allowing 
to (numerically) study QCD from its first principles. However, the price to pay is 
the appearance of diverse discretisation artifacts like breakdown of the Lorentz in- 
variance, a necessity to work with the Euclidean formulation of the theory and, in 
practice, at finite volume. We discuss in details the methods allowing to handle most 
of the artifacts. Lattice QCD has been successfully used in phenomenology (mass of 
the charm quark, B- and D-mesons physics, generalised parton distributions, QCD 
at finite temperature and its phase diagram, etc.). But it can also be used to study 
the fundamental parameters (like coupling constant) and properties of the theory 
itself. This is the main goal of the present dissertation. We concentrated our effort 
on the study of the two-point correlators of the pure Yang - Mills theory in Landau 
gauge, namely the gluon and the ghost propagators. We are particularly interested 
in determining the Aqcd parameter - the fundamental scale of the pure Yang-Mills 
theory. It is extracted by means of perturbative predictions available up to NNNLO. 
The related topic is the influence of non-perturbative effects that shows up as appea- 
rance of power-corrections to the low-momentum behaviour of the Green functions. 
We shall see that these corrections are quite important up to energies of the order of 
10 GeV. 

Another question that we address is the infrared behaviour of the Green func- 
tions (in Landau gauge), at momenta of order and below Aqcd- At low energy the 
power-law dependence of some Green functions changes considerably, and this is 
probably related to confinement. The knowledge of the infrared behaviour of the 
ghost and gluon propagators in Landau gauge is very important, because many 
confinement scenarii (for example the Gribov-Zwanziger scenario) give predictions 
for their momentum dependence at very low energies. Ab initio simulations on the 
lattice is a quasi-unique method for testing these predictions and the only way to 
challenge the underlying models for confinement. 

We try to clarify the laws that govern the infrared gluondynamics in order to un- 



derstand the radical nature of the changes of the infrared behaviour of some Green 
functions. Many questions arise : the Gribov ambiguity, the validity of different non- 
perturbative relations (like Schwinger - Dyson equations, Slavnov - Taylor identities 
and renormalisation group equations) at low momenta, self-consistency of the lat- 
tice approach in this domain. The lattice approach allows to check the predictions of 
analytical methods because it gives access to non-perturbative correlators. Our main 
goal thus is to use lattice Green functions as a non-perturbative input for different 
analytical relations. This allows to control the approximations that are done within 
the traditional truncation methods for the non-perturbative relations. Such a mixed 
numerical-analytical analysis of the complete Schwinger - Dyson equation for the 
ghost propagator provided us with an interesting alternative to the widely spread 
claim that the gluon dressing function behaves like the inverse squared ghost dres- 
sing function, a claim which is at odds with lattice data. According to our analysis 
the Landau gauge gluon propagator is finite and non-zero at vanishing momentum, 
and the power-law behaviour of the ghost propagator is the same as in the free case. 
However, as we shall see, some puzzles remain unsolved. 



Chapitre 1 

Continuum and lattice QCD 



In this chapter we recall very briefly the most fundamental ideas of QCD (defi- 
nitions, symmetries, covariant gauge fixing, Gribov ambiguity). As we are inter- 
ested in a non-perturbative calculation of different correlators, we also introduce 
diverse non-perturbative relations between Green functions, namely Schwinger- 
Dyson equations, Slavnov-Taylor identities and exact renormalisation group equa- 
tions. After this we shall discuss the lattice formulation of the pure Yang-Mills theo- 
ries, in particular the procedure of the Landau gauge fixing on the lattice. 

1.1 General features of QCD 

1.1.1 Definitions and symmetries 

Nowadays there is no doubt that Quantum Chromodynamics (QCD) is the theory 
of the strong interaction. The fundamental principle of this theory is local gauge 
invariance. This principle, together with general principles of locality, Lorentz in- 
variance and renormalisability, imposes important constraints on the form of the 
Lagrangian. The simplest form in Euclidean four-dimensional space reads 

£qcd = -\^ a]lv + E WDrf-mito (1.1) 

xp=u,d,s,c,b,t 

with (go is the bare coupling) 

F; v = d F A a v - d v A a F + g Q f ahc A h v A c v , a = l..N 2 c - 1 (1.2) 
D }l = d }t -ig Q t a A a ir (1.3) 

This Lagrangian is invariant under gauge transformations of the fields 

A v {x) ^ A^\x) = u(x)A }l (x)u f (x) + i[d }t u(x)]u A (x) (1.4) 
tp(x) i — ^ xp( u \x) = u(x)xp(x), (1.5) 

where u(x) e SU(N C ) and Nq = 3 is the number of colours. 

In order to quantise QCD using the functional integration formalism one has to 
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integrate over the quark and the gauge bosons fields. The Grassmannian integral on 
the quark fields is Gaussian, that is why we discuss only the integration on the gauge 
boson fields A. The fields A^ u ' and A in the equation 1)1 .4[) are related by a gauge 
transformation, and thus they are physically equivalent. So, in order to quantise a 
gauge theory one performs an integration over gauge transformation equivalence 
classes - the orbits of the gauge fields. This is the Faddeev - Popov procedure. The 
integration on the orbits is done by choosing a representative element on every orbit, 
i.e. fixing the gauge with some relation 

f[A] = 0. (1.6) 

The condition / should define the orbit of the field A in a unique way. Then the 
generating functional for Green functions reads 

Z\j,WM = j [VAVipVlp] A f [A] 5 (/ [A]) e $ ^Qcd+J d*x(A«j*+wf+lpa;) ^ (1 7) 

where all loop integrals are understood to be regularised. We denote the ultraviolet 
cut-off ar x , go = g(a^ 1 ). The Faddeev-Popov determinant Ay[y4] which appears in 
this formula is defined by means of invariant integration 

A f [A] J Vu(x)Y[S (f \a {u \x)^) = 1 (1-8) 
yielding in the general case 1 



6f 



= E det' 1 1 J (L9) 

t:f[A^]= °* 

Choosing the Landau gauge condition 

f[A] : dpAp = (1.10) 
and supposing for the moment that it fixes the gauge in a unique way, one obtains 

A L andau[-4] = det (A + igod^) = J [VCVC] r/WKl^Of). (1.11) 

The spurious anticommuting fields c and c belonging to the adjoint representation 
of the gauge group are called Faddeev-Popov ghosts, and 

M%{x,y) = (A + igod^A^S^ix-y) (1.12) 

is the Faddeev-Popov operator. The corresponding formula for the generating func- 
tional can be easily generalised by choosing for the gauge condition 

f[A] : BpAp = a(x), a(x) e su(N c ). (1.13) 



^^when the condition lll.6t does not fix the gauge in a unique way. 
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In this case Ay remains the same as ()l.llj> , and one can integrate on a(x) with some 
Gaussian weight having a dispersion £. This gives for the generating functional 



Z\i,Q, co, a, tr] = J [VAVtpVtpVcVc] e~ f d i xC eff [A,^,c,c]+^ i (L14) 
C eff [A,Tp f y,c,c] =£qcd- (3 ^' )2 -c a (x)(5 ab A + ig f abc A c tl d ll )c\x) (1.15) 
L = y ti 4 x + a)t/> + ipo; + crc + cc) . (1-16) 



The choice £ = corresponds to the Landau gauge. The gauge fixing term in l)1.14|) 
can be expressed as a result of Gaussian integration on an auxiliary field B a (x). This 
gives another form of the Lagrangian C e ff : 

£brst = Cqcd ~ ~ (B a ) 2 + B a d }l A a }l + c a (5 ab A - ig f abc d }i A c fi )c b . (1.17) 

The QCD Lagrangian written in this form is invariant under BRST transformations. 
If A is a constant infinitesimal Grassmann number these transformations take the 
form 

5A a u = XD a u c c c 

r r 

Sip = ig M a xp 

bc a = -\g^f abc c b c c 
5c a = XB a 

5B a = 0. (1.18) 

The virtue of the BRST transformation is its global nature. This simplifies a lot the 
derivation of the Slavnov-Taylor identities (direct consequence of the gauge inva- 
riance). We discuss this question below. 



1.1.2 The Gribov ambiguity 

A serious theoretical difficulty pointed out by Gribov (T]| arises when performing 
the quantisation of a non-Abelian gauge theory (in covariant gauge) in the case of 
large field magnitudes. The reason for this is the non-uniqueness of the Landau 
gauge condition ()1.10|l . Indeed, let us find all the intersections of the gauge orbit 
with the hypersurface defined by fll.lOD . Imposing the Landau gauge conditions for 

both fields A^ and Ay? in ()1.4|l we obtain the following equation for u(x) : 
\ A }l {x) ^ Af{x) 

I d }l A }l _=Q -> dJu(x)D ¥ (x)i?(x)) = 0. (1.19) 

{ ^Af = 
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Setting at the leading order 

u ~ I + z'a(x), w + ~ I - icL(x) oc(x) e su(N c ) (1.20) 

we obtain the following equation for a : 

dp (a^a + i[^,a]) =0 — ► fyD^a = 0. (1.21) 

But 9 ?; D^ is nothing else but the Faddeev-Popov operator in the covariant gauge. 
Thus, any non-trivial zero mode of the Faddeev-Popov operator generates an in- 
tersection point of the gauge orbit with the hypersurface dl.lOI I. If this point is not 
unique, we speak about the so called Gribov copy. All these secondary gauge confi- 
gurations correspond to the same physical field A }U and thus they must be removed 
from the functional integration measure in the partition function. 




FlG. 1.1 — The Gribov region. 

A solution to this problem is to supplement the initial gauge-fixing condition 
dl.lOI I with some additional requirement. Gribov's solution consists in restricting 
the integration measure (for the gluonic field) in l)1.14j> to the domain where l)1.21[) 
has a unique solution, see Fig. 11.11 This domain is called the Gribov region, and 
its boundary (where the Faddeev-Popov determinant vanishes) is called the Gribov 
horizon. It has been argued that some of topological solutions like instantons lie on 
this boundary Q. 

Inside the Gribov region all the eigenvalues of the Faddeev-Popov operator are 
strictly positive 2 . Hence one can realise the Gribov quantisation by using the Mini- 
mal Landau gauge. In this gauge one integrates on the fields satisfying the ordinary 
Landau gauge and belonging to the set of local minima of the integral 

J d*x (aJ(z))\ 

As we shall see, this ensures that all the proper values of the Faddeev-Popov ope- 
rator are positive. The Minimal Landau gauge will be discussed in details in the 
subsection 11.2.21 

2 we recall that in the Euclidean formulation the Faddeev - Popov operator is hermitian. 
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Nowadays it is known that the Gribov quantisation prescription is not exact - 
there are secondary solutions to ()1.19|> for some fileds inside the Gribov region. The 
domain free of them is located inside the Gribov region, and it is called the funda- 
mental modular region [3], [4]. This means that the correct quantisation prescription 



A 



det FP = 



FlG. 1.2 — Gribov region and the fundamental modular region. 

consists in restricting the integration in ()1.14|l to the fundamental modular region A 
instead of the Gribov region G, see figure 11.21 However it is argued in Q that the 
expectation values calculated by integration over the Gribov region or the funda- 
mental modular region are equal. So the Gribov quantisation prescription would 
become in fact exact. We discuss this question in details in the section l2~4l 



1.1.3 Schwinger-Dyson equations 

The Schwinger-Dyson(SD) equations is a specific class of non-perturbative equa- 
tions relating Green functions and vertices. They can be easily derived in the func- 
tional integration formalism (for a review see (6|). Let Z[J] denote a normalised 
(Z[0] = 1) generating functional l)1.14jl for the Green functions, and W[J] = logZ[/] 
- the one for connected Green functions. Then the effective action, which is a gene- 
rating functional for one-particle irreducible (1-PI) vertex functions, is obtained by 
the Legendre transformation 

Tm = *m h - Wll] , fi = *m, a.*) 

where (pi denotes the generic field (A^,c fl , . . .) and ; ; is the corresponding source. 
Then, introducing the action S = J d^xC e tt and using the quantisation prescrip- 
tion based on the integration over the fundamental modular region, the Schwinger- 
Dyson equations are obtained from the observation that 

/ [ii^,]^- s+H ) - / [n^iy s+H (-^+h) = / [n^^- s+H 

A 1 A 1 dA t 1 

(1.23) 

Using the Zwanziger's argument |5[ (quoted above) on the equivalence of integra- 
tions over A and G we restrict the integration domain to the Gribov region. It follows 



1-t 
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from this that the integral on the boundary vanishes because the Faddeev-Popov de- 
terminant present in the integration measure is equal to zero on 9G (by definition). 
This allows to write the whole set of the Schwinger-Dyson equations for the Green 
functions in a compact form 



6<Pi 



S 

L*/J 



+ h)Z[}]=0. 



(1.24) 



We use a generic relation between two smooth functions f(x) and w(x) 



f 



( d 
\dx 



M x ) — p w ( x ) 



f 



dzu(x) d 
dx dx 



1, 



(1.25) 



that can be applied to ()1.24|> and yield the equations for the functional W generating 
the connected Green functions 



8S[cp] \8_ + SW 



1 + ji = 0. 



Scpi [8j 5] 

Finally, performing a Legendre transformation l|1.22)> . we have 



(1.26) 



8S\£ 
Sep, 



8 2 W 5 

H 

Y 8 8 5<p 



8(pi 



(1.27) 



corresponding to the full set of Schwinger-Dyson equations for proper functions. 



As an example, we derive explicitly the SD equation for the full ghost propagator. 
Varying S with respect to the antighost field c a we obtain 

Varying the last relation with respect to o~ b (y) one obtains 

{w(x) Ch{y) ) =^ (4) (^-y) = ((d p D^(x))c\y)). (1.29) 
Denoting the full ghost propagator as 

F^ ah (x,y) = (c a {x)c h {y)Y (1.30) 

we obtain ( A(x,y) = 5^\x — y)5 is the Laplace operator) 

S^(x-y) = A(x r z)F^(z,y)+ig d { fl x) (A }l (x)c(x)c(y)) . (1.31) 

Note that this equation can be obtained in a simpler way only using the definition 
of the Faddeev-Popov operator. Indeed, the ghost correlator in the background of a 



given gauge field configuration = A^t a is given by 

F i ( conf (A^y) = X^^y)" 1 . (1.32) 

The subscript means here that the equation is valid for a given gauge configuration. 
Thus one obviously has 

5\x - y) = Mjf^^g^A^), (1.33) 

where a summation on z is understood. Using the explicit formula ()1.12|> for AApp 
we get 

S(x-y) = A(x,z)F^(z f y)+ig d^ (a h (x)F^ (A,x,y)) , (1.34) 

valid for any gauge field configuration A. Performing the functional integration on 
A one gets the mean value on gauge configurations 

8{x - y) = A(x,z) (F^fry)) + igotf ] (A F (x)F^(A,x,y)) . (1.35) 
Using ^Fj^ nJf (z,y)y = F^ 2 \z,y) we find the equation dl.31l . The translational inva- 

(x) -\(v) 

riance of the Green functions allows to replace o„ by — V . Performing the Fourier 
transform on [x — y), we have finally 

1 = -p 2 F (2 V) - igoP? (MWSL^p)) ■ (1-36) 

This derivation elucidates the trivial dependence of this equation on the functional 
integral weight with which we calculate the average (•) on the gauge fields A. The 
form (11.341 1 allows an explicit discussion of the Gribov copies dependence of the 
solutions. We address this question below. 

Performing the Legendre transformation for the three-point function in (11.311 ) 
and introducing the ghost-gluon vertex 

Tf c (-q,k;q- k) = gof abc (iq v i)T v i v (-q,k;q - k) (1.37) 

and the full gluon propagator G^J ab (p), we write the Schwinger-Dyson equation for 
the ghost propagator in Fourier space 

( Fi2) Tab (k) = ~ Khkl " Sof aCd f bf J -^-J^mq v ')fM-q,kq ~ k) (ik pi ) (q - k), 

(1.38) 

given in a diagrammatic form at Figure IT31 This equation is much simpler than the 
one for the gluon propagator, because the last involves complete three- and four- 
gluon vertices (cf. Figure ll.4|) . Another virtue of ()1.38|> is that its form is explicitly 
independent of the choice of the integration domain in the functional integral (11.231 ), 
because the equality (11.341) holds valid for individual gauge configurations. 



1U 
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o- 



-1 




FlG . 1.3 — Schwinger-Dyson equation for the ghost propagator in a pure gauge theory, diagrammatically. 







G 



FlG. 1.4 — Schwinger-Dyson equation for the gluon propagator in a pure gauge theory, diagrammatically. 



1.1.4 Slavnov-Taylor identities 



The Slavnov-Taylor identities (|7J,|8|) are Ward identities in the case of a non- 
Abelian gauge theory. These identities follow directy from the gauge symmetry. We 
write the generating functional in the case of a general gauge fixing condition (11.61) 



Z\j] = J [DA] detM FP e 



(1.39) 



and use the fact that the functional integration measure is invariant under the spe- 
cific gauge transformations 



y\abJo 
LJ y, fc 



Sf a [A] =M%e h . 



(1.40) 



The second equation is just a general definition of the Faddeev-Popov operator. The 
integration measure times the Faddeev-Popov determinant is invariant with respect 
to ()1.40|l , and hence 



Z[j] = j [DA] detMppe 
Choosing a particular form for e 



(1.41) 



M 



FP 



be 



(1.42) 



we obtain a functional relation for the generating functional Z 



r 



s_ 



be 



5_ 



M 



-l 

FP 



x,y 



'Sj 



m = o. 



(1.43) 



In principle, one can use the relation l)1.25J> in order to obtain an equation for the 
functional W[j]. However, this derivation would lead to a very cumbersome expres- 
sion. In fact it is much easier to use the Slavnov-Taylor relations obtained within 
the BRST formalism. The main idea of the derivation remains the same because 
the BRST transformation is just a specific form of the gauge transformation J1.40I I, 
M pp (x, y) being the propagator of the Faddeev-Popov ghost in the background field 
A. One has using lll.17ll.18J> and the BRST invariance of the generating functional 

m 



d 4 x J [VAVcVc}e-S dixC ?ff 



■ D a "c 



ab Jj 



2 J 



abCgbgC 



0. (1.44) 



This equation (J1.44JI allows to obtain the Slavnov-Taylor identities for the Green 
functions by differentiating with respect to the sources j, cr, a. Writing ()1.44Jl in terms 
of generating functionals W and T, one obtains (see |9J, [7J, [8J for details) the gene- 
ral form of Slavnov-Taylor identities between propagators and proper vertices. The 
relation that we shall use in the following relates the three-gluon vertex T\„v if' °l> r ) 
to the ghost-antighost-gluon vertex, and involves the complete ghost and gluon pro- 
pagators. It reads (Q/Ell) 



r 2 G( 2 )(r 2 ) 



V-Vv r Afi (r 7 p;q)- 



(1.45) 



-lAlu) T\ v {c\,p;r 



Some remarks regarding the non-perturbative validity of the Slavnov-Taylor iden- 
tities are in order. The above derivation is invalid when M.pp is singular (see (II .401 ), 
l)1.42J> and ()1.43Jl ), i.e. for gauge fields lying on the Gribov horizon. However, this 
transformation is well defined inside the Gribov horizon. Note also that the general 
form of the Slavnov-Taylor identities does not depend on the choice of the integra- 
tion domain inside the Gribov horizon (A or G). Another argument in favour of 
non-perturbative validity of the Slavnov-Taylor identities may be given within the 
stochastic quantisation formalism 0. 



1.1.5 Renormalisation group equations 

For the sake of completeness we present here another set of non-perturbative 
relations between the correlators, namely the exact renormalisation group equations 
(or ERGE, see II Oil for a review). Those are flow equations describing the variation 
of the effective action with the infrared (or ultraviolet) cut-off. The infrared cut-off 



1U 
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is introduced by adding a special term to the action 

AS /c = \YLj ^<Pi(P)MKp 2 )<p t (-p), (1-46) 

where the momentum cut-off function for the field (pi satisfies 

r a,(* /P 2)->o p>k 

\Uk,pi)^k* P <k. {lA/) 

The role of AS^ is to suppress quantum fluctuations with momenta below k. Then 
one may show that the partition function satisfies the equation 

^ = -I I wr^^ijmk)- (1 - 48) 

The problem that arises within the formalism of the RG equations is the loss of the 
gauge invariance caused by the cut-off term. However, this problem can be solved 
by considering a modified set of Slavnov-Taylor identities 1111 . One may express the 
equation ()1.48|l in terms of the generating functional for proper vertices, that lead to 
an infinite system of partial differential equations relating different Green functions 
and the cut-off function ljl.47j> . We shall review some of the results for solutions of 
the truncated system of such equations in the section |4~T1 



1.2 Lattice QCD 

1.2.1 Lattice QCD partition function 

A fully non-perturbative study from the first principles of QCD phenomenon re- 
quires a direct calculation of the functional integral of the type ()1.14|l . These integrals 
can be approximately evaluated by means of lattice simulations. Another interest in 
the lattice regularisation is that it preserves the gauge invariance. The inverse lattice 
spacing a 1 plays the role of the ultraviolet cut-off, and we recover the continuum 
limit theory by sending a to zero. 

In what follows we discuss only pure Yang-Mills theories. In practice, when 
doing a lattice simulation, one considers a theory in a finite volume V = L 4 with 
(most often) periodical boundary conditions ; and generates some (quite large) num- 
ber M of gauge field configurations {Q} distributed according to the probability 
measure 

dp. [A] = e~ J" ^Yang-Mills W [£>^] . ( L49 ) 

Then one can calculate a Monte-Carlo approximation for any operator O 

/-i M 
dpO(A v ) ^ - Y,0(Q). (1.50) 
i=l 

Let us now discuss the measure dp[A] in the discrete case. Gauge fields are defi- 



ned on the links of the Euclidean lattice (cf. Figure 1X31) , and the fundamental lattice 
variable for the gauge field is the link variable 

U(x,x + ejl ) = U }l {x) = e^Mx+W) e SU(N C ), (1.51) 
where e« is a vector in direction ]i, || = a. For small values of the lattice spacing a 




FlG. 1.5 — Gauge fields A% are defined on the links of the lattice (||e f( || = a), and the fermion fields ip(x) are 
defined on its sites. 

we can extract the field Au using an approximate formula 

2ia g0 ' =A(* + f) + W . (1.52) 

We use this definition of the gluon field in what follows. Then one can define an 
elementary gauge invariant variable - a plaquette 

U(p) = U^U^x + ae^Ulix + aeMix). (1.53) 

Using this variables we can write the simplest action that converges to the pure 
Yang-Mills' action in the continuum limit ( 1121 . see 1131 , 111411 for a review) : 

^[^W]=ffiE(l-^TrU(p)), (1.54) 

where the lattice bare coupling is defined by 

= (1.55) 

So 

Thus, the partition function in the lattice formulation reads 

ZlatM = J [VU }t (x)}e- s *M. (1.56) 
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Of course there exist an infinite number of lattice actions giving ()1.14|l in the conti- 
nuum limit a — > 0. They should all give compatible results when the lattice spacing 
is small enough and the renormalisation group scaling laws are verified. In practice, 
one generates the gauge configurations {Q}, z = 1 . . . M according to the probabi- 
lity measure [Ptyx^e-V^.Inall our simulations we have used the Wilson action 
J1.54I I. Usually, the Metropolis or Heatbath algorithms are used for the Monte-Carlo 
generation. For a detailed review on this topic see |13|.|14|. 

The important question that has not been discussed up to now is the removal 
of the ultraviolet divergences of the theory in the continuum limit. The regularised 
lattice partition function contains one parameter - the bare coupling constant go (or, 
equivalently, the lattice spacing a). When calculating some physical quantity on the 
lattice, say, the string tension o~, we obtain it as a function of a. If the ultraviolet cut- 
off a -1 is large enough, the calculated quantity obeys the scaling law, and one can 
compare the result for a (a) of the lattice simulation to the known experimental data 
(7 ex p, and thus determine the corresponding physical value of the lattice spacing a by 
solving the equation a (a) = o~ ex p. In the unquenched case one has to calculate several 
physical quantities because more non-fixed parameters are involved (the masses 
of quarks). When all free parameters of the lattice theory are fixed in the scaling 
region with some experimental inputs, all further calculations 3 are automatically 
renormalised, and do not contain any divergences. Moreover, all the calculations 
are now the predictions of the theory. Thus, the only limitations of the numerical 
method a priori are discretisation errors and the errors on the experimental inputs. 
We discuss in details the systematic error of lattice simulations in the section l23l In 
practice one is also often limited by the computer power. 

Let us say some words about the strong coupling limit of (|1.56|) . corresponding 
to the perturbative expansion in the f> parameter. The ordinary perturbation expan- 
sions in gauge theories (in powers of go) are at most asymptotic, but power series 
of (11.5611 in /3 is proven to have a finite radius of convergence II15II . Many interes- 
ting results, like the existence of the mass gap and the area law for the Wilson loop, 
have been analytically proven within this approach. However, it has been argued 
that the region of the strong coupling is analytically disconnected from the weak- 
coupling domain, corresponding to the continuum limit of the theory 111 61 . Lattice 
gauge theories are believed to possess an essential singularity at some finite value 
trough/ corresponding to an infinite order phase transition. This phenomenon is well 
known in statistical physics (roughening transition). Thus it is impossible to know 
whether the strong-coupling predictions are applicable in the physically interesting 
weak-coupling regime. In practice, one has to perform numerical simulations with 

$ > trough- 

1.2.2 Fixing the Minimal Landau gauge on the lattice 

The continuum gauge fixing condition dl.6l l is modified by the discretization, so 
one works with its lattice version /l(U^) = 0. A gauge configuration Uq generated 
during the simulation process does not satisfy the Landau gauge condition. One has 
to perform a gauge transformation u(x) on it in order to move the field configura- 



3 of experimentaly observable quantities 
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tion along its gauge orbit up to an intersection with the surface fi(Uu) = 0. But 
there is no need to have an explicit form of u(x). Instead we perform an iterative 

minimisation process that converges to the gauge fixed configuration U^K Let us 
first illustrate it on the example of the Landau gauge in the continuum limit. For 
every gauge field A we calculate the functional 

F A [u(x)} = -Tr J d 4 x (a^ 1 \x)A { f' '(x)) , u(x) e SU(N C ). (1.57) 

Expanding it up to the second order around some element uq(x) we have (u = e x UQ, 
X e su(N c )) 

F A [u] = F A [u ] -2 J (fixTr (xa^j" o) ) + j d 4 xTr (x f M T ? [^ ( " o) ] x) + . . ., 

(1.58) 

where the matrix Aijsp in the quadratic term defines the Faddeev-Popov operator. 
Obviously, if uq is a local minimum of (11.571 1 then we have a double condition 



M-pF 



is positively defined. 



(1.59) 



Hence, the minimisation of the functional 111 .571 ) allows not only to fix the Landau 
gauge, but also to obtain a gauge configuration inside the Gribov horizon (because 
all the eigenvalues of .Mfp are positive). On the lattice, the discretized functional 
\1.57\ reads 

F u [u(x)} = -— 9*eTr^w(x)lI F (x)tt + (x + e F ). 



(1.60) 



X,\i 



Then at a local minimum uq we have a discretized Landau gauge fixing condition 



E 

/' 



A 



(«o) 



x + f)-A ji 



that we write in a compact form V^A^ = 0. Indeed, if u(x) 
co(x) is the element of the algebra su(N c ) then 



("o) 



X ~2 



6_ 

5u 



F(u) 



"0 



d 
ds 



F(u(s,x)) = 0, 



(1.61) 

Uo(x)e scv ^ where 
(1.62) 



s=0 



and hence 

1 

~ V 



9=leTr^ ([o;(x + e F ) - a>(x)] U H (x)) = Va>(x) G su(N c ). (1.63) 

X,]l 



Thus 



Yj Uh(x) - ilj(x) - U pi (x - e }l ) + ilj(x - e ¥ ) 
V 



= 



(1.64) 
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and at leading order in the lattice spacing a dl.52ll we obtain (I1.61I) . The second deri- 
vative (the equivalent of the second-order term in l)1.58|l ) can be written as 



d 2 1 / 

j^F(u(s,x)) = -- [to^pAy 



(1.65) 



where 



A ¥ = 1 (-cv(x)Up\x) + U { ; i \x)cv(x + e F ) + co(x + e^U^ (x) - (x)co(x) 

(1.66) 

This defines a quadratic form 

(to,Mfp[U}(vj = -^eTr'£ j ^[co 2 (x + e ll )-2co(x + e F )co(x)+co 2 {x)} uf°\x] 



2V 



x,n 



Tr£ ( {w{x + ey) - co{x)f (u^\x) + U^ ] \x 



r(«o)+/ 



(1.67) 



The operator M^p[U] is the lattice version of the Faddeev-Popov operator. It reads 



Mf P [U]cv) a (x) = ^S [ Gf(x) 



co (x + e») — co (x) ) — (x <-> x — e 



where 



G;\x) = ~Tr({t\t b }(u fi (x) + U^ 

a h { x + y) = Tr ( fC (*W - "£(*))) • 



(1.68) 

(1.69) 
(1.70) 



For the Minimal Landau gauge fixing in our numerical simulation we have used 
the Overrelaxation algorithm |17| with co = 1.72. We stop the iteration process of 
the minimising algorithm when the following triple condition is fullfield : 



V(N} - 1) 



X,}i 



v„<'">) (v^<" w » 



V(N} - 1) 



Ag(*!)-Ag(fe) 



< ®max x \d p A«\ = 10 



< e ju = io- 



< &sA n = io- 7 . 
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(1.71) 



Ag(*0+Ag(fe) 

where iM\x) is the matrix of the gauge transformation u(x) at the iteration step n, 



and the charge 




(1.72) 



must be independent of t in Landau gauge when periodical boundary conditions for 
the gauge field are used. The choice of numerical values for the stopping parameters 
is discussed in |18J. 



Chapitre 2 

Lattice Green functions 



In the functional integral formalism Green functions are defined as mean va- 
lues of products of fields according to the functional measure, giving as a result the 
vacuum expectation values of these products. Often one is interested in the Green 
functions in Fourier space. Here we define the Fourier transformation on the lattice. 

If a function f(x) is defined on the sites of the four-dimensional lattice with per- 
iodical boundary conditions then 

~ In 

f(p)=a 4 LK x ) e ~ tP ' X Vv = ^n v , n,=0,l,...,L-l 

/(*) = ^E^' x /(p)- 

V 

In the case of the variables defined on the links of the lattice one should change the 
above formulae by x« — ► x« + % . In the infinite-volume limit V — > oo we have 




We recall in the first part of this chapter the main definitions regarding the lattice 
two- and three-point Green functions in Landau gauge and the Momentum sub- 
straction renormalisation scheme. Next we describe the details of the numerical cal- 
culation of the ghost propagator and discuss different sources of errors of the lattice 
approach. The last part is devoted to the Gribov ambiguity on the lattice and the 
influence of Gribov copies on the Green functions. 



2.1 Green functions in Landau gauge 

The gluon propagator in Landau gauge may be parametrised at all values of 
momenta as 

d$«\ V ,- V ) - (Aj(-p)AJ(p)) = S* (V-^) G (2) (P 2 )- (2.3) 



This implies that the scalar factor (p 2 ) may be extracted according to 

G(2)(p2) = 3(4-1) ^^"^"^ p2 * ° (z4) 

completed with 

G^(0) = 4(N , 1 _ l) E4f a (0^)^ p 2 = 0. (2.5) 

The difference in normalisations at zero ()2.5|) and finite momenta l)2.4[) is due to an 
additional degree of freedom related to global gauge transformations on a periodic 
finite lattice (p 2 = <-> x„ ~ + L). 

The ghost propagator is parametrised in the common way : 

F (2)ab( p/ _ p) = (?(-p)f(p)} = 5 ab F®(p 2 ). (2.6) 

It is not defined at p 2 = because in this case the Faddeev-Popov operator is strictly 
equal to zero and thus it is not invertible. 

The renormalisation scheme that is widely used in order to renormalise the lat- 
tice Green functions is the so-called Momentum substraction scheme (MOM). The 
virtue of this scheme is its non-perturbative definition. The renormalisation constants 
are defined by setting the corresponding Green functions to their tree values at some 
renormalisation point p 2 . In the case of the two-point Green function the renorma- 
lisation constant Z 3 (p 2 ) of the gauge field or the one of the ghost field (denoted 
Z3(p 2 )) is unambiguously defined as : 

Z 3 (p 2 ) =p 2 G^(p 2 ) (2.7) 
Z 3 (p 2 ) =p 2 F (2) (p 2 ). (2.8) 

The coupling constant has also to be renormalised to complete the renormalisation 
of a pure Yang-Mills theory. It can be defined non-perturbatively by an amputa- 
tion of a three-point Green-functions from its external propagators. But this requires 
to fix the kinematic configuration of the three-point Green-function at the norma- 
lisation point. On the lattice one usually uses either a fully symmetric kinematic 
configuration (denoted MOM) or a zero point (ZP) kinematic configuration with 

one vanishing external momentum (denoted generically MOM), see Figure l2~Tl 



Three-gluon vertex : symmetric case 

There are only two independent tensors in Landau gauge in the case of the sym- 
metric three-gluon Green function 1191 : 

r £,W/' 3 (Pi' P2, Ps) = <W 2 (Pi - n)n + <Ws(P2 - ps) Fl + <W(P3 - Vi)n ( 2 -9) 

T [2] /_ _ _ x _ (Pl-P2) m (p2-P3h 1 (p3-Plh 2 
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Then the three-gluon Green function in MOM scheme (p\ — v\ = f \ = V- 1 ) can be 
parametrised as 



+ 



(2.11) 



For a non-perturbative MOM definition of the renormalised coupling gR one need 
to extract the scalar function G^ s ^ m (fi 2 ), proportional to the coupling go at the tree 
order. This is done by the following projection : 



i fabc . ^ ^ ^ 



(2.12) 



Three-gluon vertex : asymmetric case 

The three-gluon Green function with one vanishing external propagator (T20I. IT9I ) 
can be parametrised as 

(2.13) 

and thus 

G (3)as ym(p2) = ^ _Z_ VppG y ( p /0/ _p). (2.14) 



Ghost-gluon vertex 



Similar parametrisation may be written in the case of the ghost-ghost-gluon 
Green function (cf. Figure I2.2|> . But in this case one obtains two different renor- 
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FlG. 2.1 - Definitions of the MOM scheme, t> 2 = fi 2 (left) and MOM, q 2 = r 2 = k 2 = pi 2 (right) 



malisation schemes MOM c and MOM c o corresponding to the zero-point kinematic 

configuration with vanishing momentum of, respectively, the gluon and the ente- 
rs) 

ring ghost. We denote by Gj> a generic scalar function extracted from a three-point 
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FlG. 2.2 — Kinematic configurations oftheMOM c , MOM c q schemes, p 2 = ]>i 2 



function in a particular kinematic configuration K. Then the gauge coupling at the 
renormalisation scale }i 2 is defined by 

in the case of three-gluon vertices, where the choice of p, determines the renormali- 
sation scheme. For of ghost-ghost-gluon vertices the coupling is defined by 

where GjjP {jp\, p\, p^) is the scalar factor of the ghost-gluon three point function 



*4f((Pi) c (P2)c(p3)y. In fact, only the coupling in the MOM c scheme can be defined 
non-perturbatively for the reasons explained in the following section. 



2.2 Numerical calculation of the ghost Green functions 
in Landau gauge 

2.2.1 Lattice implementation of the Faddeev-Popov operator 

In order to calculate numerically ED the ghost propagator 




(2.17) 



one uses the lattice definition (|1.68|> of the Faddeev-Popov operator. Most lattice 
implementations of the Faddeev-Popov operator use its explicit component form. 
But as we have seen the action of the Faddeev-Popov operator on a vector co can 
also be written as a lattice divergence 

Mf P [U]cv = (2.18) 

where A' is defined by (I1.66I) . This form allows a very efficient lattice implementa- 
tion which is based on the fast routines coding the group multiplication law : 



# a>j n , co ou t are the ghost fields. 

# Uu{x) is the gauge configuration. 
type SU(N C ) U fl (x),dU,W,W + ,W- 
typesu(N c ) co in (x),co out {x) 

for all x do 
dU = 0. 
W = cv in (x) 
do fi = 1 ... 4 

W+ = a>;„(x + e /( ) 

W_ = w in (x-e }l ) 

dU = dU + U fl (x - e }l ) * W + W * U }l {x) 

-U li (x)*W+-W-*U p (x-e li ) 

end do 

w out {x) =dU- dU f - -^-MdU - dU f ) 
end do 



2.2.2 Numerical inversion of the Faddeev-Popov operator 

We invert the Faddeev-Popov operator by solving the equation 

J^Mf P [U] ab (x,y)f ] b (y)=S a (x), (2.19) 

y,b 

for some source S a (x) using an appropriate algorithm (for a review of algorithms see 
(22|). The operator _Mpp[li] has zero-modes, that is why the inversion can only be 
done in the vector subspace orthogonal to the kernel of the operator. The trivial 
zero-modes are constant fields. If we neglect Gribov copies then the Faddeev-Popov 
operator has no other zero-modes, and thus the non-zero Fourier modes form a basis 
ofX^ : 

V (y)= £c/>*, WrjeK'-. (2.20) 



The inversion in one Fourier mode 

Choosing the source for inversion in the form 

SJ(x) = <J 8 V>* 



(2.21) 



and taking the scalar product of the inverse M l pp[U]~ 1 Sp with the source one obtains 

(SJ | Mfp-ty) ) = E ( (^"X) e " <P ' ( *" y) (Z22) 



= yf( 2 )(p) (2.23) 

after averaging over the gauge field configurations. This method requires one matrix 
inversion for each momentum p. It is suitable only when one is interested in a few 
values of the ghost propagator. 



The inversion in all Fourier modes 

One can calculate the ghost propagator for all momenta p doing only one matrix 
inversion noticing that 

s *» = h + h £ e ~ iHx ' y) (Z24) 

V P7^0 



and choosing for the source 



S a Q (x)=S ab [5 Xi0 - v ). (2.25) 



1 



xy 



The Fourier transform of M 1 Sq, averaged over the gauge configurations, yields : 

= E^ f(2, w4E^ If(2) (^-y) 

x v x,y 

= F^-\p)-5{p)Y J ^ 2 \x), (2.26) 

X 

where we have used the translational invariance of the ghost propagator. Therefore, 
with this choice of the source, only one matrix inversion followed by one Fourier 
transformation of the solution is required to get the ghost propagator for all values 
of the lattice momenta. 

Because in the case of the source ()2.25|) one inverts in all modes at the same time, 
some statistical accuracy is lost. However, it turns out to be sufficient for our pur- 
poses. 

There is one important technical point that should be mentioned. During the in- 
version process it is mandatory to check, whatever the choice of sources, that roun- 
ding errors during the inversion do not destroy the condition that the solution is 
still orthogonal to the kernel of the Faddeev - Popov operator : 

E (Mfp'h) (x) = (2.27) 



Indeed, if the zero-mode component of the solution grows beyond some threshold 
during the inversion of the Faddeev-Popov operator on a gauge configuration, then 
this component starts to increase exponentially, and a sizeable bias is produced in 
other components as well. We have observed this phenomenon occasionally, about 
one gauge configuration every few hundreds, when using the componentwise im- 
plementation of the lattice Faddeev-Popov operator based on J1.68I I. However, we 
have never observed sizeable deviations from (|2.27|> using the efficient implementa- 
tion of the Faddeev-Popov operator exposed in the subsection 12.2.11 



2.2.3 Calculation of the ghost-gluon vertex 

In order to calculate the ghost-gluon three-point function in MOM c and MOM c o 
renormalisation schemes one has to calculate the corresponding ghost two-point 
function : 

MPiWVifciVi)) latti f E ^MiM^m-^ip^Ps) (2.28) 

v ' conf. i 

It is quite easy to calculate this Green function in the MOM c kinematic configura- 
tion, because in this case pj_ = —p, p$ = p and (A4pp[Uj] _1 ) (— p, p) is just a ghost 
propagator in the background gluon field defined by 

The situation changes when considering the kinematic configurations like MOM c o, 
when the momentum of the entering (or of the outgoing) ghost is set to zero. In this 
case the inversion of the Faddeev-Popov operator has to be performed with the 
source ()2.21|) . In other words we try to solve the equation 

Mf P [U] xyVyz = e l ^ x - z \ (2.29) 

The vector rj y — yY^ z t]y Z e ip ' z is the solution of 

M\%[U] xy rj y = fr*, (2.30) 

and 

V y = ^LVy^ P ' Z = ^c(y)c(zyv z = c{yfc{p). (2.31) 

z v z 

Doing a summation on y we obtain a Fourier transform of the field c(y) : 

^LVy = c(0Mp)- (2-32) 

y 

But the last equation expresses the orthogonality condition ()2.27|> . Thus we find 
c(0)c(p) = in this case. That means that ()2.28|> is also zero, and the vertex func- 
tion cannot be directly extracted (on the lattice) in the kinematic configuration with 

vanishing ghost momentum MOM c o, but only in the MOM c scheme. 



2.3 Errors of the calculation 

There are three main sources of errors when calculating Green functions on the 
lattice : the statistical errors of the Monte-Carlo method, the systematic bias coming 
from the space-time discretisation and finally the error due to the gauge fixing (in- 
fluence of lattice Gribov copies). The last is discussed in the following section. 



2.3.1 Estimating the statistical error 

The gauge-field configurations produced via the Monte-Carlo generation pro- 
cess (see [13J for a review) are not completely decorrelated. However, the residual 
correlations may be neglected. Nevertheless, all data points (as function of momen- 
tum p) of a Green function are calculated on the same set of gauge configurations 
{Q}, and in this sense they are not independent. This problem arises when calcula- 
ting quantities involving functions of mean values, like the coupling constant in a 
MOM scheme. In order to take in account the bias induced by this correlation one 
uses a special method (called Jackknif e 11231 , 11241 ) of computation of the error. Gene- 
rally speaking this method is a standard bootstrap method (of the estimation of the 
variance in the case of a non-Gaussian distribution) based on a resampling with re- 
placement from the original sample. We start with a Monte-Carlo sample of size M. 
Our purpose is to calculate the error on the estimation of the mean of this sample. 
We divide it into [M/ m] groups of m elements : 



0(C 1 ) f ... f O(C n 



0(C m+ i),. . .,0(C2 m ) 



...,0(C 



M 



(2.33) 



Next one defines the partial averages 



E£i^(Q)-E 



a k 



(k+l)m 
i=km 



0{Q 



k = l...M 



M-m 

and finally obtains the following expression for the error : 



M 



m 



(2.34) 



^jackknif e 

(O) 



\ 



M - 1 
M 



/ 



M 

E' 

k=l 



Lfc=l a k 
M 



(2.35) 
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This analytical expression differs from the standard formula for the dispersion of 
the mean value by an additional factor ~ M. 



2.3.2 Handling the discretisation errors 

Because of discretisation of the space-time, lattice theory 1 looses the rotational 
symmetry SO (4) inherited from the Lorenz invariance in Minkowski space. This 
symmetry is replaced by a discret isometry group H4 = S4 x P4 (semiproduct of the 



we suppose that the lattice is hypercubic. 



permutation and reflection groups) having 4! • 2 4 = 384 elements. A generic scalar 
function G(p) extracted form Green functions is thus invariant along the orbit 0(p) 
generated by the action of the group H4 on the components of the lattice momentum 
p = x (tii, n 2> n 3> n 4)- It m ^y be proven in the theory of group invariants that each 
orbit is characterised by four group invariants 

pM=fl B J>J, n = 2,4,6,8 (2.36) 

V- 

One may average on the gauge orbit of the H4 group in order to increase the statis- 
tics : 

A5 l(p P!,pW, pM/P H ) = ^^ E% (2.37) 

where card O(p) is the number of elements in the orbit 0(p). The resulting average is 
a function of the four invariants p'"'. But in the continuum limit any scalar function 
is a function of the rotational invariant p^. We will explain how it is possible to 
remove the dependence on the three other invariants with the example of the free 
lattice gluon propagator : 

G (p) = -X^, where ^ = ?sin(^). (2.38) 
l^fi V}i a \ a / 

If all the components of the lattice momentum verify the condition apu <C 1, then 
yj^ p 2 ~ p 2 — Y2_a 2 p^, and thus one has up to terms of order ~ a 4 

+ 0(fl 4 )j. (2.39) 

So, taking the continuum limit a — > is equivalent to taking the limit p^ —>■ 0. We 
apply this idea to ()2.37|l . Making the reasonable hypothesis of regularity near the 
continuum limit, we expand : 

G L (p [2] ,p [i] ,p [6] ,p [8] ) « G L (p[ 2 U,0,0) + P [4] py(p [2 l0,0,0) + ■ ■ • (2.40) 

and Gi(p^, 0,0,0) is nothing but the scalar factor in the continuum in a finite vo- 
lume, up to lattice artifacts which do not break 0(4) invariance. When several orbits 
exist with the same p 2 , we can remove an important part of the hypercubic artifacts 
by extrapolating the lattice data towards Gi(p^, 0,0,0) using a linear regression 
with respect to pi 4 ' (the other invariants are of higher order in a). 

There are vectors p that have only one orbit. In order to include them in the data 
analysis one should interpolate the slopes obtained in the extrapolation of ()2.40|> . 
This can be done either numerically or by assuming a functional dependence of 
the slope with respect to p 2 . In principle the p 2 dependence of the slope may be 
calculated using lattice perturbation theory. A dimensional analysis suggests the 
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FlG. 2.3 — Example of the extrapolation in p^for the ghost scalar factor F (p) in the case of the SU(3) gauge 
group for the lattice volume 32 4 and at f> = 6.4. Red round dots correspond to the bare data, and blue squares 
to the extrapolated one. In practice we do not consider the moments above ap = S. 



form ~ : c } 2 - We have used the function 



if)' 



dG 



(p [2 U0,0) = -i-2 (c 1+ c 2 pW) (2.41) 



dpW ' ' ' ' ( pP] ) 



with two fit parameters in order to fit the slopes. The validity of the exposed me- 
thod is qualitatively checked by the smoothness of the resulting curve. At Figure l2~3 
we present an example of removing the hypercubic artifacts. We see that the me- 
thod works very well, even at large values of ap. In practice we do not consider the 
momenta above ap = 5. 



2.4 Gribov ambiguity and lattice Green functions 

We have discussed in the previous sections different uncertainties introduced 
by the discretisation of space-time. Another bias comes from the gauge fixing (see 
section ll.l.2|l . As we have already mentioned, the Minimal Landau gauge fixing 
quantisation is equivalent to realising the Gribov quantisation prescription. Thus 
lattice Green functions are calculated within this prescription. However, the gauge 
is not fixed in a unique way. The Gribov ambiguity on the lattice shows up by the 
non-uniqueness of the minimum (11.591) of the functional (11.601 . Indeed, the functio- 
nal (11.601) has a form similar to the energy of a spin glass which is known to have 



exponentially-many metastable states (cf. Figure EU). All these lattice Gribov copies 
are situated inside the Gribov horizon. On a finite lattice one can (in principle) fix the 
gauge in a unique way However this is very expensive in terms of computer time. 



energy of a spin glass 




spin configurations 



FlG. 2.4 — Landscape of minima for a generic spin glass. Circles mark the configurations corresponding to 
metastable states. Their number grow like exp (const ■ N) as function of the number of spins N. 

Nevertheless one has to understand the influence of the choice of the minimum of 
the functional dl.601) on the mean values yielding Green functions. For this purpose 
we studied the landscape of minima of the gauge-fixing functional 1181 . In the follo- 
wing subsection we define a specific probability to find a Gribov copy, as function 
of the physical momentum. In the next-to-the-following subsection we discuss the 
influence of these copies on the Green functions, and the Zwanziger's conjecture 
on the equivalence of the integration over the Gribov region and the fundamental 
modular region in the infinite volume limit. 

2.4.1 The landscape of minima of the gauge-fixing functional 

Let us consider the landscape of the functional Ffj. One of its characteristics is 
the distribution of values at minima F m { n of F(j. We know that for small magnitudes 
of the gauge field all the link matrices U^(x) (they play the role of couplings bet- 
ween the "spin" variables) are close to the identity matrix, and thus the minimum is 
unique. Fheir number increases when the bare lattice coupling jS decreases, because 
the typical magnitude of the phase of LL(x) grows in this case and thus link ma- 
trices move farther from the identity matrix. Fhe number of minima also increases 
with the number of links (at fixed /3) because in this case there are more degrees of 
freedom in the system. 

We can define a probability to find a secondary minimum, as a function of the 
jS parameter. For each orbit we fix the gauge Ngf times, each gauge fixing starts 
after a (periodic) random gauge transformation of the initial field configuration. We 
thus obtain a distribution of minima F m { n . This distribution gives us the number of 
minima N(F l ) as a function of the value of F l = Fi- . The relative frequency of a 
minimum F l is defined by 

.. N(P) 
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FlG. 2.5 — Probability (averaged over gauge orbits) to find a secondary minimum as a function of f> at 
volumes V = 8 4 , 10 4 , 12 4 , 16 4 . Solid line represents a fit with an empirical formula \2Ab\ . The vertical dashed 
lines delimit the window of each fit. 



where £j N(F 1 ) = Ngf- Then the weighted mean number of copies per value of F m i n 
is given by 

N = Y j cv 1 N(F'). (2.43) 

i 

This allows us to define a probability to find a secondary minimum when fixing the 
Minimal Landau gauge for a given gauge field configuration : 

N 

= 1 - g^y. (2-44) 

If one finds the same value of F m i n for all Ngf tries then this probability is zero. On 
the contrary if all F; are different then pi con f is close to one. 

Having the probability to find a secondary minimum when fixing the gauge for 
one particular configuration we can calculate the Monte-Carlo average (•) on gauge 
orbits, i.e. on "spin couplings" U^{x). We obtain finally the overall probability to 
find a secondary minimum during a numerical simulation at given /3 : 
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We have performed simulations in the case of SU(2) latice gauge theory at vo- 
lumes V = {8 4 , 10 4 , 12 4 , 16 4 } for f> varying from 1.4 to 2.9. For each value of jS we ge- 
nerated 100 independent Monte-Carlo gauge configurations, and we fixed the gauge 
Nqp = 100 times for every configuration. Between each gauge fixing a random 
gauge transformation of the initial gauge configuration was performed, and the mi- 
nimising algorithm stops when the triple condition (|L7TJ) is satisfied. Examples of 
the resulting probability to find a secondary minimum are presented at Figure l2~5l 

As expected, the probability is small when /3 is large, and it is close to one when /3 
is small. The dispersion was calculated using the Jackknife method (discussed in the 
subsection l2.3.1l l. The physical meaning of this dispersion is the following : when the 
error is small, all gauge configurations have a similar number of secondary minima. 
On the contrary, this dispersion is large if there are some exceptional gauge configu- 
rations having a different number of copies. At small /3 almost all gauge configura- 
tions have many secondary minima, that is why the dispersion of the probability is 
small. At large /3 almost all gauge configurations have a unique minimum, but some 
of them can have copies. This may considerably increase the dispersion of the proba- 
bility. The appearance of exceptional gauge configuration possessing a large density 
of close-to-zero eigenvalues of the Faddeev-Popov operator has been recently re- 
ported II25I . Probably these fields are related with those having a lot of secondary 
minima at large /3's, and this correlation deserves a separate study. 

We can fit our data (Figure I2.5|> for the probability ()2.45jl with an empirical for- 
mula 

p C) = < 2 ' 46 > 

in order to define a characteristic coupling f> c when the probability to find a copy de- 
creases considerably. One can define /3 C as corresponding to the semi-heights of the 
probability function P(/3). At this value of /3 an equally probable secondary attractor 
of the functional Fjj appears. The fit has been performed for the points between da- 
shed lines at Figure 12.51 and the results for the fit parameters are given ibidem and 
in Table EU We see that /3 C depends on the volume of the lattice. Let us check whe- 
ther these values correspond to some physical scale. According to works |26"1.||2"7| 
one has the following expression for the string tension a for /3 > 2.3 : 

[era 2 ] 06) ~e h ^ Uo ' hP (2.47) 



with c = 4.38(9) and d = 1.66(4) 2 . Using this formula, we define a characteristic 
scale corresponding to the critical values /3 C from Figure 1231 : 

A c = fl(jSc) • L 

in the string tension units, L is the length of the lattice. In the last column of the Table 
I2.1l we summarise the results. We see that for the values of /3 C in the scaling regime, 
when the formula ()2.47j> is applicable (/3 > 2.3), we obtain compatible values for the 
physical length A c . 



in this cited formula we kept the original convention for the RG /5-function, namely j6g 
ll/3N c , ft = 17/3N£. They are different from the one we use in the following. 



L 


ft 


1 T 

^/ndf 


A c , in units oil/ ^/a 


8 


2.221(14) 


0.27 


3.85(31) 


10 


2.342(6) 


0.28 


3.20(21) 


12 


2.44(1) 


0.22 


2.78(20) 


16 


2.541(5) 


1.92 


2.68(17) 



TAB. 2.1 — TTze characteristic length defining the appearance of secondary minima. The errors for A c include 
errors for d and c parameters, and the fitted error for f> c 



This suggests that in the case of SIT (2) gauge theory lattice Gribov copies appear 
when the physical size of the lattice exceeds a critical value around 2.75/ \fcr. At first 
approximation A c is scale invariant, but a slight dependence in the lattice spacing 
remains. 

In principle, the parameter B c can be calculated with good precision. One should 
do it in the case of SU(3) gauge group, because the dependence of the lattice spacing 
on the bare coupling is softer, and the scaling of the theory has been better studied 
than in the case of the SU(2) theory. 

A natural question that arises after the study of the distribution of F m [ n is whe- 
ther the gauge configurations having the same value of F m i n are equivalent, i.e. they 
differ only by a global gauge transformation. This can be checked by calculating the 
two-point gluonic correlation function on the gauge configuration. Indeed, accor- 
ding to the lattice definition of the gauge field that we used dl.52l l, 

G fconf(* -y) « Tr [( u A x ) - *■£(*)) ■ (^(y) - u t(y))] • ( 2 - 48 ) 

Applying a global gauge transformation U^(x) — > VU^(x)V^ we see that the gluon 
propagator remains unchanged. This is also the case of the ghost propagator scalar 
function. We have checked numerically that the values of the gluon and the ghost 
propagators in Fourier space are the same for gauge configurations having the same 
F^. Thus we conclude that gauge configurations having the same F m i n are in fact 
equivalent |18J. 



2.4.2 Lattice Green functions and the Gribov ambiguity 

As we have mentioned in the subsection ll.1.21 there is a conjecture |5f saying that 
in the infinite volume limit the expectation values calculated by integration over the 
Gribov region or the fundamental modular region become equal. Let us briefly recall 
the argument given in |5j. Let denote a field lying on the Gribov horizon : 



f a F B F = o 

\ Mfp(B)coq = 0, cvq(x) ^ const, coq(x) e su(N c ). 



(2.49) 
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We write Fg(£,o;) for the functional (11.571 ), with u(x) = exp tco(x). On the boundary 
one has 

F' B (0,cv ) =0 
F£(0,o, )=0 

F7(0,o;o) = |/d 4 x [d F (o;2)] 2 cxy 
So for small variations of t one has for the functional dl.571) 

F B (t) = F B (0) + ^'(0)* 3 + ^"(0)t\ (2.51) 

and the last expression is minimised for t = t = —3F'J'(0) / Fg"(0). The value that is 
achieved at this secondary minimum inside the Gribov horizon is 

F B (t) = F B (0) - - [ B J 3 . (2.52) 

Using the estimations (12.501) and the fact that the third order derivative appears in 
an even power, one sees that the secondary minimum is lower than the one on the 
boundary, and the difference between them decreases with the volume V. The corres- 
ponding gauge configuration can be written as 

B ¥ {x,l) = B }l (x,0)+t [D ¥ {B)co Q ] (x) = B H (x f 0)-3^j^ [D }t (B)cv ] (x). (2.53) 

When the lattice volume V is large, the integration (over G or A) on the gauge confi- 
gurations is dominated by a small shell near the boundary, because the dimension 
of the configuration space is large (2dV(N^ — 1)). According to the above argument 
the configurations near the boundaries 9G and 9 A draw together. So, having (12.531) 
for a typical gauge field configuration, it is natural to suppose that in the infinite 
volume limit the average calculated by integration over the domains G or A be- 
come equal. However at a finite lattice this is clearly not the case, that is why it 
is very important to know what is the influence of lattice Gribov copies on the 
Green functions. This question that has already been considered by different au- 
thors (EHMZQ/EDMSIMSZMSHI/IZSI ). To check the dependence of Green functions 
on the procedure of the choice of the minimum we adopted the same strategy as in 
above citations : for every of the 100 gauge configurations used to compute Green 
functions the gauge was fixed 100 times (a periodic random gauge transformation is 
done after each gauge fixing). The Monte-Carlo average was computed with respect 
to the "first copy" (fc) found by the minimisation algorithm and the "best copy" (be), 
having the smallest value of F m i n . We have calculated the gluon and the ghost pro- 
pagators, and also the three-gluon Green functions in symmetric and asymmetric 
kinematic configurations. The simulations have been performed in the case of the 
SU(2) group on lattices of volumes 8 4 and 16 4 for /3 = 2.1,2.2,2.3. According to the 
results of the subsection l2.4.1t at these values of /3 we are sure to have lattice Gribov 
copies. We conclude Ifl8ll that no systematic effect could be found for gluonic two- 
and three-point Green functions, the Monte-Carlo average values in the cases of (fc) 
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0.211 
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0.045 


2.1 
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-16 


0.145 
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0.078 
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16 
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-16 


0.023 
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2.3 
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1 




0.086 
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0.024 


2.3 


16 


4 


„[«] = 


-16 


0.114 


0.034 



TAB. 2.2 — Volume dependence of the ghost propagators JTSf, P;< = w« 



]6 


L 


(fmin) { u] 


^(fmin){tJ} 


2.2 


8 


-0.8236 


0.003744 




10 


-0.8262 


0.002367 




12 


-0.8272 


0.001377 




16 


-0.8279 


0.000802 


2.4 


8 


-0.8642 


0.005270 




12 


-0.8669 


0.002739 




12 


-0.8686 


0.001849 




16 


-0.8702 


0.001003 



TAB. 2.3 — Volume dependence of the Monte-Carlo+gauge orbit mean value at minima F m i n and the disper- 
sion of this mean. 

and (be) being compatible within the statistical errors. However the ghost propaga- 
tor is quite sensitive to the choice of the minimum - in the case of (be) the infrared 
divergence is lessened. This dependence has been found to decrease slowly with 
the volume [33 J. The results of the subsection 12.4.11 indicate that the convergence 
can happen only beyond the critical volume A 4 . To check this we compare the fc/bc 
values of the ghost propagator, at one physical value of the momentum, for the orbit 
n 1 = 1 on a 8 4 lattice and the orbit n 2 = 4, n' 4 ' = 16 on the 16 4 lattice 3 at the same f> 
(see Table l2~2)> . It happens indeed that the decrease is observed only at /3 = 2.1 and 
2.2, in accordance with Table 12.11 However, these values of /3 are not in the scaling 
regime, and thus a study on larger lattices would be welcome. It is not surprising 
that the ghost propagator depends on the be/ fc choice : the (be) corresponds to the 
fields further from the Gribov horizon where the Faddeev-Popov operator has a zero 
mode, whence the inverse Faddeev-Popov operator (ghost propagator) is expected 
to be smaller as observed. The correlation between the bc/fc choice and the gluon 
propagator is not so direct. Another quantity is obviously strongly correlated to the 
bc/fc choice : the value of F m [ n . We tested the volume dependence of the Monte- 
Carlo+gauge orbit mean value of the quantity F m [ n (see Tab l2.3l) . According to the 
argument given in |5|, all minima become degenerate in the infinite volume limit, 
and closer to the absolute minimum (in the fundamental modular region). We see 
from the Table l2~3l that their average value and dispersion decrease with the volume 



3 Remember that the momentum in physical units is equal to Inn/ '(La) in our notations. 



— VJ11UUV aiALU'l^UltJ til LU lCiLLl\_C U1CC1L 1U1LLUU110 



at fixed f>, in agreement with 0. 

We finish this section with a brief summary of our results Ifl8ll regarding the Gribov 
ambiguity on the lattice : 

1. Lattice Gribov copies appear and their number grows very fast when the phy- 
sical size of the lattice exceeds some critical value 2.75/ y/cr in the case of 
the SU(2) theory). This result is fairly independent of the lattice spacing. 

2. The configurations lying on the same gauge orbit and having the same F m { n 
are equivalent, up to a global gauge transformation, and yield the same Green 
functions. Those corresponding to minima of Fjj with different values of F m [ n 
differ by a non-trivial gauge transformation, and thus they are not equivalent. 

3. We confirm the result (OHI/ESI/EDI/EII/ESI/ESI/EI) that the divergence of 
the ghost propagator is lessened when choosing the "best copy" (correspon- 
ding to the choice of the gauge configuration having the smallest value of Fjj). 
We also showed that gluonic Green functions calculated in the "first copy" and 
"best copy" schemes are compatible within the statistical error, no systematic 
effect was found (with periodic gauge tranformations). 

4. We found that the influence of Gribov copies on the ghost propagator de- 
creases with the volume when the physical lattice size is larger than the critical 
length discussed above. We also show that the quantity F m i n decreases when 
the volume increases. These two points are in agreement with the Zwanzi- 
gers's argument [5j on the equality of the averages over the Gribov region and 
the fundamental modular region. 



Chapitre 3 

The ultraviolet behaviour of Green 
functions 



The Lagrangian of the Pure Yang-Mills theory in a four-dimensional space-time 
does not contain any dimensional parameters susceptible to fix an energy scale for 
dimensionless quantities. However, the spectrum of the corresponding quantum 
theory contains massive states (glueballs). As a matter of fact, the quantum theory 
possesses a finite energy scale called Aqcd/ which is generated by the quantisation 
process followed by the renormalisation. All dimensionful physical quantities are 
expressed as multiples of powers of this scale, and thus it should be a renormalisa- 
tion group invariant : 



A 



QCD 



o, 

(3.1) 



where /3 (g(}i 2 )) is the renormalisation group beta function and ji is the renormali- 
sation scale. The solution of the above equation reads 



Aqcd [H,g(l* 




(3.2) 



where g\ is an arbitrary integration constant. Aqcd is a renormalisation scheme- 
dependent quantity, although it is a renormalisation group invariant within one 
particular scheme. So it is not a real physical quantity. Still, its value is important 
for estimating the lowest bound of the domain of validity of perturbation theory. 
Knowing several first coefficients of the /3-function we find from the equation ()3.2|l : 



A, 



QCD 



}iexp 



1 



2/3o \gj g 2 (^ 



+ 



log 



g\ 



+ 0(g 2 ). (3.3) 



We see that there is an essential singularity when g 2 (}i 2 ) — > 0, and thus a perturba- 
tive calculation of related quantities (for example, the string tension yfa = c^Aqcd) 
is impossible. In the following sections we describe the method of calculation of 
Aqcd from lattice Green functions in Landau gauge. We start with a review of the 
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purely perturbative results for Green functions. The momentum range available on 
the lattice is situated at rather low energies where the non-perturbative power cor- 
rections are not negligible. The section I5~2l is devoted to the estimation of the domi- 
nant power corrections. At the end of this chapter we present the results of analysis 
of our lattice data. 



3.1 Aqcd and perturbative expressions for Green func- 
tions 

Different scalar factors of Green functions depend on the Aqcd parameter dis- 
cussed in the previous subsection. These scalar factors can be calculated non-perturbatively 
in lattice simulations, and one can extract Aqcd by fitting the lattice data in the ultra- 
violet domain to the corresponding perturbative formulae. Here we make a review 
of available perturbative Landau gauge calculations for the ghost and gluon propa- 
gators in the MOM schemes. 

If {pi,g\, ji 2 ) is a renormalised proper vertex in Landau gauge, then the cor- 
responding proper bare vertex function is independent of the renormalisation point 
ji. This fact is reflected by the Callan-Symanzik equation for the renormalised func- 
tion : 

(aEJ? + ^ ( gRi}i2) ) T g -^ (W 2 ))) ^\rugW) = o 0.4) 

where 7 (^r(^ 2 )) is the anomalous dimension. In the Momentum subtraction schemes, 
the renormalisation conditions are defined by setting some of the two- and three- 
point functions to their tree-level values at the renormalisation point. Then (13.41 1 
simplifies to 

d\n(Z 3MOM (p 2 = }i 1 ,ar x ) 
lim — — 2 = 73,mom(#mom) (3.5) 

in the case of two-point Green functions, where Z^(ji 2 ) is defined in ()2.7|> , a 1 stands 
for the ultraviolet regularisation and 73,mom(^mom) is the anomalous dimension. 
A similar expression can be written for the ghost propagator renormalisation factor 
Z3. As we have already seen in the section IZll there is an infinite number of MOM 
schemes differing by kinematic configurations at the substraction point. We limit 
ourselves to the configurations defined by the subtraction of the transverse part of 

the three-gluon vertex (MOM) and that of the ghost-gluon vertex with vanishing 

gluon momentum (MOM c ) and vanishing incoming ghost momentum (MOM c o), 
discussed in the section IZT1 

Both anomalous dimensions for ghost and gluon propagators have been recently 



computed ( 11341 , 113510 in the MS scheme. The result at four-loop order reads 
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(3.6) 



where h = g 2 / (An) 2 . In order to obtain the coefficients of the anomalous dimensions 
1)3 .5|l in a MOM scheme one has to express the above expressions in terms of the 
corresponding coupling. We use the results of the article Il34l where the three-loop 
perturbative substraction of all the three-vertices appearing in the QCD Lagrangian 
for kinematic configurations with one vanishing momentum are given. In Landau 
gauge and in the pure Yang-Mills case one has the following relations between the 
couplings in different MOM schemes and : 
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(3.7) 



Thus, inverting ()3.7J> and substituting in H3.6J). we obtain the anomalous dimensions 
of the gluon and ghost propagator in the three above-mentioned renormalisation 
schemes. In a MOM scheme, the equations (|3.6|) may be integrated as functions of h 
(cf.l53t 1 : 
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■"■We omit the index specifying the renormalisation scheme both for /z and Aqcd in the following 
formulae 
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+ /z 3 (/3 3 73 - /3 2 7 2 + (/3 /3 2 - /3 2 /3 2 ) ^ 
+ (-/3 2 /3 3 + 2/3 /3 1 /3 2 -/3 3 ) 7 )-l I + 



(3.8) 



where % are the expansion coefficients of the anomalous dimension in a generic 
MOM type scheme and Zq is an integration constant. The knowledge of the f>- 
function 



dh 



(3.9) 



at some order n allows to calculate the momentum dependence of h. At four-loop 
order one has 



+ 



^te + KI) 3 (- 2log3(f)+5log2( ' )+ ( 4 - 6 if) 1081 



0-i 

(3.10) 



where t = In — . The last equation together with <|3.8h allows us to write the ghost 



QCD 



and gluon propagators as functions of the momentum and Aqcd- The numerical 
coefficients for the /3-f unction in H3.91 are summarised in the Table |3~T1 : 
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MOM c0 




11 


h 


102 
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2412.16 


2952.73 


3040.48 


h 


84353.8 


101484 


100541 



TAB. 3.1 — The numerical coefficients for the ^-function for different MOM schemes 136F 



3.2 OPE for the Green functions and dominant power 
corrections 

The momentum dependence of the QCD Green functions at low energies is mo- 
dified by non-perturbative effects. These effects show up by presence of power- 
corrections to logarithmic series or, in other words, by non-zero values of corres- 
ponding condensates. For example, such a non-perturbative object as instanton has 

a weight <x exp — frrr/ giving a * leading order a power correction oc \. It is argued 

s (P ) p 

in ([37J,[38J) that non-perturbative lattice gluonic two- and three-point functions 
include such contributions up to quite large energies of around 10 GeV. For a sys- 
tematic study of Aqcd one nas to know the influence of power corrections on the 
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Green functions. 



A powerfull tool to study the dependence of Green functions on the non-pert- 
urbative condensates is the Operator Product Expansion (OPE) II39H . This method is 
applicable to the problems having a specific energy hierarchy, or two very different 
characteristic energy scales. For example, in QCD it may be applied to the study 
of the influence of some background semi-classical field configurations. We recall 
here the idea of this method on the example of a two-point correlation function of a 
generic field (p 



G(x) 



.2/ T V 2) I 

It is postulated that when x — > the product of the fields may be expanded as 



♦(§)♦(-! 



EIXMof-'to), 

w=0 i 



(3.12) 



\2n] 



where the second sum is performed on all local operators 0\ of mass dimension 
2n having the same quantum number than the l.h.s. The OPE suggests that all the 
features of the short-distance behaviour are stored in the Wilson coefficients 



x 



series m a s 



(3.13) 



that can be calculated in perturbation theory. In Fourier space they behave as 



and thus 



wf(p) ~ -j 



00 

g(p) = E E 5 ? K' lo g ( v L iy 



/ i \ 




w x 


series in a s 



n=0 



,-1 



(9 



In 



,w+l 



(3.14) 



(3.15) 



where the coefficients vf are computed in perturbation theory, and (oj 2 "') are va- 
cuum condensates. At n = 0, corresponding to the trivial basic operator I, we find 
an ordinary perturbative series for G. But other condensates may lead to the ap- 
pearance of non-perturbative power corrections. Usually this method is applied to 
gauge-invariant product of currents, and involves only gauge invariant quantities 
(for a recent review see 11401 ). However it can be extended to gauge-dependent ope- 
rators (like QCD propagators) and involve gauge-variant condensates ( 1411 . 1421 ). We 
do not discuss here the subtile question of the renormalisation of condensates and 
of calculation of their anomalous dimensions. On the lattice the MOM-type renor- 
malisation process is non-ambiguous (1 431. 114411 ). because the non-perturbative value 
for the l.h.s in (13.151) is available. This allows to define the condensates at fixed ultra- 
violet cut-off. Then one can apply a MOM renormalisation prescription on the both 

sides of ()3.15|> and thus renormalise the condensates (oj 2 "'). 

In the following paragraph we will discuss the dominant power corrections, and 
corresponding condensates, in the case of the gluon and ghost correlators. 
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3.2.1 The dominant OPE power correction for the gluon propaga- 
tor 

The basis of operators in the pure Yang-Mills case is 

J_ c a d }l A a v c a c b A a F c b A a }l A b v d F c a c a c b c a c b ... (3.16) 



At the leading order (a oc 1/p 2 power correction compared to perturbation theory) 
only underlined operators contribute l43l to the gluon propagator, because opera- 
tors with an odd number of fields cannot satisfy colour and Lorentz invariance and 
cc does not contribute because of the particular structure of the ghost-gluon vertex 
(cf. Figure l3~Tl b)). We write then for the gluon propagator : 




FlG. 3.1 — (a) Contribution of the gluon A 2 — condensate (represented as soft external gluons) to the gluon 
two-point function (b) Contribution of the ghost cc condensate (represented as soft external ghosts) to the gluon 
and ghost two-point functions. These contributions vanish because they are proportional to (^zero) momentum 
of the outgoing ghost is the ghost-gluon vertex. 



(g (2) )^(p 2 ) = (t(a;(- p )a^(p) 



• A c (0)A C (0) - \ (3.17) 

/ T/ \ao /„z\ , / T/ \aba'b' ( Jl\:a'b': V pK ' pK ' 7 , 
K v 0)jiv \V > + WjtvjiV VP ) d d u>v> 4( ^ N 2 _ fj + • • •' 

where (•) is a v.e.v with respect to the non-perturbative vacuum and : A£(0) A° p (0) : 
is a free-field normal product. In the perturb ative vacuum the v.e.v. of all the normal 
products give zero, and thus only Vq is non-vanishing. Hence 



(Vo);Up 2 ) = (g^YJP 2 )- (3-18) 



ab 

I 

fiV 

The coefficient V2 is obtained at the tree-level order from 



7> v ' 9 

and 



(g\:A c J0)A c J0)-.\g)=2 + O(a s ) (3.19) 



(g\T (K(-p)K(p) ) |g> connected = {V 2 )iit (p z )(g\ : AJ,(0X(0) : \g), (3.20) 



UVfl 
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where \g) is a soft gluon state. So, using the LSZ rule to cut the soft external gluons, 
we obtain 

M u 2 1 (^r(0)A« (-p)AJ(p)A«(0)) 



(g^(O)) " ,(g<&(0) 



which can be computed in perturbation theory (cf. Figure Ell a)). Finally, 

(3.22) 

A MOM-type renormalisation prescription may be defined non-perturbatively. This 
allows an easy renormalisation procedure for the A 2 — condensate |43(. Here we do 
not include the effects of the anomalous dimension of the A 2 operator II45II and hence 
we apply the MOM prescription by imposing the tree-level value to the Wilson co- 
efficient at the renormalisation point p 2 = }i 2 for the last equation. This allows to 
factorise the perturbative gluon propagator giving finally 

Z 3 ( M 2 ) = Z 3 ,p er t(M 2 ) fl + y f^fl^ + Oiglr 4 )) ■ 0-23) 



3.2.2 The dominant OPE power correction for the ghost propagator 



In the case of the ghost propagator the set of basic operators is the same, the 
ghost condensate cc does not contribute for the same reasons as for the gluon pro- 
pagator (cf. Figure I3~l7 b) ). Thus, applying the OPE to the ghost two-point function, 
we obtain : 

F (2) "V) = WV) + (Vl)"*" (p 2 )(: K(°)K(o) ■) +■■■ 
where, in analogy with (13.211 1, the Wilson coefficient reads 

, i , ;^>-m^,(o)t( c ^) ^',(o)) connected 

w ah = (V 2 ) S s % r = ~ S j-± '-, (3.25) 



which is equal to twice the diagram represented on the Figure I372l that describes the 
coupling of the ghost propagator to the gluon A 2 — condensate. Hence 

w° b = l -5 s % r ■ 2^ (i go f^) ^ (igofW) ^ = N c p£:\f). (3.26) 



U-l LI CI V 1U1C L Id V 1U Lll Ul \J1 CU I X U.X IA_ L1U1 ID 




FlG. 3.2 — Contribution ofthegluon A 2 — condensate (external soft gluons) to the ghost propagator 

This gives the leading non-perturbative contribution, because the first Wilson coef- 
ficient trivially gives the perturbative ghost propagator. Finally, 

F (2) "V) = ?<2f (V 2 ) (l + ^ jfl^y + O (^5,,- 4 )) (3.27) 

where all quantities are bare. Performing the MOM renormalisation we obtain for 
the renormalisation factor : 

z 3 ( F 2 ) = z^T t (F 2 ) ( 1 + ^ 4 fif! ) i ) +°(gi^ 4 )) ' ( 3 - 28 ) 

where the A 2 -condensate is renormalised as in the case of the gluon propagator. We 

see that the dominant multiplicative correction to the perturbative Z3 ;pert is identical 
to the one obtained in the previous section for the gluon propagator (I3.23I I. 

3.2.3 Constraints on the Wilson coefficients from the Slavnov-Taylor 
identity 

The gauge-dependent power corrections due to the (A 2 ) -condensate are obviously 
absent in gauge-invariant quantities. Because of this the Wilson coefficients for the 
(A 2 )-condensate in different Green functions are not independent. Some relations 
may be obtained from the Slavnov-Taylor identity ()1.45|) but their role in the MOM 
renormalisation constants is not obvious. 
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FlG. 3.3 — The (A 2 ) contribution to the ghost- gluon vertex with vanishing entering momentum. The above 
diagram is zero in Landau gauge because of the projector in the gluon propagator. 

It is interesting to know if there are any power corrections to this vertex, be- 
cause perturbation theory predicts the non-renormalisation of this vertex, i.e. it is 
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FlG. 3.4 — Non-zero dominant (A 2 ) contribution to the ghost-gluon vertex with vanishing entering momen- 
tum. This term contribute to the external ghost propagator. 

equal to 1 to all orders (00/[34J). If the (A 2 ) power corrections are present they will 
constitute the main contribution at low energies. One can directly evaluate the Wil- 
son coefficient to the ghost-gluon vertex T a ^ c ( — p, 0; p ) with vanishing entering ghost 
momentum. 

The only non-zero correction to the ghost-gluon vertex with one vanishing ghost 
momentum is the one that contributes to the external ghost propagator (see Fi- 
gure EU- But all the diagrams with the condensate interaction attached to different 
external legs are zero in Landau gauge (see Figure l3~3b . Thus the ghost-gluon vertex 
in this particular kinematic configuration does not contain the power-corrections, 

and thus the non-renormalisation theorem holds at this order. However, this it is not 
true if the external ghost momentum is not exactly zero. 



3.3 Data analysis 

We calculated, using the techniques described in the chapter |2j the ghost and 
the gluon propagators of the Landau gauge SU(3) gauge theory at different lattice 
volumes and different values of the /3 parameter (cf. Table l3~2)l [21]. The lattices mar- 





V 


a- 1 (GeV) 


^phys (fm 4 ) 


# Configurations 


-> 6.0 


16 4 


1.96 


6.73 


1000 


6.0 


24 4 


1.96 


33.17 


500 


-> 6.2 


24 4 


2.75 


8.43 


500 


-> 6.4 


32 4 


3.66 


8.85 


250 



TAB. 3.2 — Lattice setup parameters. The lattice spacings are taken from Table 3 in f46j with a physical unit 
normalised by \fa = 445 MeV. The lattices marked by the "— > " symbol have similar physical volume. 

ked by the "— symbol correspond to similar physical volume. The produced data 
allow us to study the propagators in the momentum range [« 2GeV, ps 6.5GeV] . 

This section is organised in the following way. In the first subsection we present 
the fits of the ghost and the gluon propagators separately, and compare the fitted 
values for Aqcd- Thus we test the self-consistency of the method. Non-perturbative 
effects are quite important in the energy interval accessible to us. This is why ano- 
ther motivation is to study the asymptoticity of the perturbative series. The lat- 
ter is done by comparing the results in different renormalisation schemes (R — 

MS, MOM, MOM c , MOM c o) and at different orders (from two to four loops). In the 
second subsection we use the analytical result of the previous section namely that 
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the dominant non-perturbative effects are the same for the ghost and for the gluon 
propagators, and hence the ratio of the gluon and the ghost dressing functions is bet- 
ter described by perturbation theory at low energies. We shall see that lattice data 
support this claim. 



3.3.1 Fitting the gluon and the ghost propagators 

We extracted Aqcd from the dressing functions of our lattice propagators by 
fitting them to the formula l|3.8|l (with h given by l|3.10|> ) in different MOM renor- 
malisation schemes. There are two parameters of the fit - the wanted Aqcd anCi the 
integration constant Zq. An example of such a fit if presented at Figure 1531 The ob- 
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FlG. 3.5 — Extrapolated lattice data at f> = 6.4 for G(p) (left) and F(p) (right). The solid line is the fit at 
four-loop order in the MS scheme. The vertical dotted lines delimit the window of each fit. 

tained value of A# in a scheme .R is converted to the MS scheme using the exact 2 
asymptotic formulae 



A 



MS 



0.346 A- 



MOM 



A ms = 0-429A~ 



(3.30) 



A 



MS 



0.428 A. 



MOM, 



c0 



The results in MS are given in Tables I3.3I3.4I331 The errors include the statistical 
error, extrapolation errors and the bias due to the choice of the fit window. We 
see from these tables that at a given order and in a given renormalisation scheme 



2 A relation between the values of Aqqd in different schemes A and B reads 



A 

A^ =6XP 



2/?o Ui(^ 2 ) s|(f 2 ) 



(3.29) 



If gg = g\ f 1 + £ + . . . J then the asymptotic freedom gives A B = A A e J Po . Thus the exact conver- 
sion coefficient is given by an one-loop calculation 1 47 1 . 





A (3) 


A (3) 
MS, ?/zosf 


MS,#Zhom 


MS,ghost 


6.0 
6.2 
6.4 


519(6)^ 

r— / y \ -1-1 / 

509(6) 27 
476(7)^ 


551(12)5* 

550 (8) 
549(7)1| 


441(4)1« 

^ /■> r- / S- \ -1-11 

435(6)1^ 
410(4)1| 


461(10)5* 
465(8)5* 
468(7) 5§ 



TAB. 3.3 — Three-loop and four-loop physical values of A^ in MeV extracted from fits in the MS scheme. 





A (3) 

MS,ghion 


A (3) 
MS,<fhost 


MS,#Zwon 


MS,ghost 


6.0 
6.2 
6.4 


324(2)5 

320(2)1? 4 

312(1)5 5 


322(8)52 
326(5)1^ 

331(4)1| 


320(4)5 


331(8)1 2 1 

353(9)5 8 7 



TAB. 3.4 — Three-loop physical values of Aj^ in MeV converted from fits in the MOM scheme. 





A (3) 

MS,%\uon 


A (3) 
MS, ghost 


Ail 

MS,#/ 


ion 


Ai 4 ! 


6.0 


345(3)5 


369(9)5 








6.2 


341(2)5 


364(8)1}! 


344(4); 


f9 
-6 


357(10)1? 6 


6.4 


323(2)1^ 


354(8)l2g 


332(2); 


1-14 

-30 


351(8)1| 



TAB. 3.5 — Three-loop physical values of Aj^ in MeV converted from fits in the MOM c scheme. 



the values obtained from the gluon and ghost propagators are consistent within the 
error bars, and are quite independent of the ultraviolet cut-off. The results from a 
direct fit in the MS scheme (Table l3~3)l confirm the old claim that we are still far from 
asymtoticity in the considered momentum interval in this scheme II48II . Our analysis 
suggests that the perturbative series become asymptotic at the NNLO in the case of 

MOM and MOM c renormalisation schemes. However, the property of asymptoticity 
is only approximate at considered momenta. To see this one can use the perturbative 
expression (analogue to (13.31) ) for A# in terms of the coupling to the order four 



21nAW-m„2 1 ^, /3p/3 2 - /3 2 ftgft - 2 ftftfr + f 
r -W< -jfa--fi ln W l R)-—fi—hR ^4 h R' 

(3.31) 

A (n+1) 

and plot the ratio (Figure l3~6b of the consecutive orders R {n) . There is a qualitative 

A R 

agreement between the ratios presented at Figure I3T61 and our results (see Tables 5-10 
in EH). The influence of truncation, responsible for the differences between different 
orders and renormalisation schemes, is mostly due to the large value of the effective 
coupling at considered energies |2TJ. In fact, as shown in l37l . the real value of the 
coupling constant may be smaller, because of the power correction discussed in the 
section I5~2l Indeed, according to the OPE analysis the effective coupling constant is 
modified by a factor 

a s -> a s ^1 + const • ^r^j ■ ( 3 -32) 
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FlG. 3.6 B j^y- for n = 2 (dashed lines) and n = 3 (solid lines), for the gluon propagator in the MS 

scheme (a), MOM scheme (c) and MOM c scheme (e), and for the ghost propagator in the MS scheme (b), 
MOM scheme (d) and MOM c scheme (f). 

According to the results of the section 13.21 one can eliminate the dominant power 
correction by considering the ratio of the propagators. In this case one expects a 
better behaviour of perturbative series at low momenta. We discuss this strategy in 
the following subsection. 



3.3.2 Fit of the ratio 

Given that at the leading order the non-perturbative power corrections factorise 
J3.23I I, (13.281 1 and are identical in the case of the ghost and gluon propagators, we can 
fit the ratio 

Z^ 2 ,A R/ (A 2 )) ggt^AR) 

Z 3 (qZ,A R ,(A*)) Z 3fPat (q2,A R )' { " ; 



iJt\_UUlL ^J.^J — UQIQ CIA IClly jlj 



to the ratio of perturbative formulae in scheme R given by d3.8D , and then convert Ar 
to Ajyjg using (|3.30|) . It is interesting to notice that non-perturbative corrections cancel 
out in this ratio even in the unquenched case with njr ^ flavours of dynamical 
quarks. The AQCD _ P arame ter extracted from this ratio is free from non-perturbative 
power corrections up to contributions related to the operators of dimension four. 
In Table l3~6l the best-fit parameters for the three schemes are presented and we plot 

in Figure 13.71 the lattice data and the MOM best-fit curve for the ratio ()3.33|) . In 




k (GeV) 

FIG. 3.7 - Plot (in the MOM scheme) of the f^T for the best fit parameter A m = 269(5) MeV. 



scheme 


MS 


r7n.d.f 


MS 


ptVn.d.f 


. 4 loops 
MS 


X 2 /n.d.f 


MOM 


324(6) 


0.33 


269(5) 


0.34 


282(6) 


0.34 


MOM c 


351(6) 


0.33 


273(5) 


0.34 


291(6) 


0.33 


MOMcO 


385(7) 


0.33 


281(5) 


0.34 


298(6) 


0.33 



TAB. 3.6 — The best-fitted values of Aj^g for the three considered renormalisation schemes. As discussed in 
the text, MOM c seems to be the one shoiving the best asymptotic behaviour. 

Figure 13.91 we show the evolution of the fitted parameter Aj^g when changing the 
order of perturbation theory used in the fitting formula. One can conclude from 

Figure l3~9l that the MOM scheme at three loops gives the most stable result for A|^g. 
It can also be seen from the ratio of four to three loops contributions (see Figure l3~T0l) 
for the perturbative expansion of In Z3, 

In Z 3 = r \nh R + r { h\ , (3.34) 
2=1 

where the coefficients r ; are to be computed from those in equations (|3.6H3.10|> using 
the Table |3~T1 The same is done for lnZ3. 
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k . (GeV) 

min v ' 



FlG. 3.8 — The determination of the optimal window fit (from 3 GeV to k max a < tz/2) results from the 
search for some "plateau" o/Axjg when one changes the low bound of the fit window. Fits are done in the 

MOM scheme. 



According to our analysis, and in agreement with the result of the separate fit, 
three loops seems to be the optimal order for the asymptoticity 3 . Indeed, the va- 
lues of Aj^g for the three considered renormalisation schemes practically match each 
other at three loops (see Figure l3~9)l . The approximate value 

A m = 269(5) + J (3.35) 

could be presented as the result for the fits of the ratio of dressing functions to per- 
turb ative formulae. 

The results of the previous subsection and (21.1 suggest that our present syste- 
matic uncertainty may be underestimated (narrowness of the momentum interval, 
truncation of the perturbative series, etc.), that is why we prefer simply to quote 
Ajyjg « 270 MeV for future reference. This value is pretty smaller than the value of 
« 330 MeV obtained by independent fits of dressing functions (see Tables l3~4l3.5l) . In 
light of our OPE analysis and previous results ||44|. this argues in favour of presence 
of low-order non-perturbative corrections to the ghost and gluon propagators in the 
momentum range [2 GeV, 6 GeV] . 



3.3.3 Comparing the results 

We showed that perturbation theory is quite successfull in describing (up to 
NNNLO) lattice propagators in the momentum range [2 GeV, 6 GeV], yielding com- 
patible values of Aqcd- The separate fit of the ghost and gluon propagators, and 
the fit of their ratio favours the existence of oc X power corrections and validates 
the OPE approach in the case of ghost and gluon propagators. The difference with 



Note that the the asymptoticity property is better verified in the case of the ratio, see Figure lS^I 
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FlG. 3.9 — Evolution of the parameter Aj^, extracted from fits of the ratio \3.33\ and propagators alone 
(rhombus and star markers, extracted from Tables 13.413.51 ET ) to perturbative formulae, as function of the 
order of perturbation theory. The solid line corresponds to the value 13.351 . Only statistical error is quoted. 
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FlG. 3.10 — (b) Ratio of four-loop to three-loop contributions (and of three-loop to two-loops for the sake 

of comparison) for the perturbative expansion of log Z3 and log Z3 (in MOM) in 13.341 plotted versus the 
momenta inside our fitting window. 



previous approaches is that we have not introduced any additional fit parameter, 
and have only used perturbation theory. Our method can also be used to calculate 
Aqcd from propagators alone in the unquenched case. In principle, it can be used 
to estimate the value of the (A 2 ) condensate. 

The main limitation of the application of perturbation theory to lattice Green 
functions in the accessible energy interval is the lack of asymptoticity and the trun- 
cation of the series. In fact, even the conversion formula J3.29I I is not exact at consi- 
dered momenta. We estimate the accuracy of our results at around 10%. It can be im- 
proved by performing the simulations at /3 = 6.6, 6.8 on the lattices of sizes 48 4 ,64 4 , 
respectively. The choice of parameters is motivated by the necessity to have the same 
physical volume of the lattice in order to control the finite-size effects. 
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Chapitre 4 



The infrared behaviour of Green 
functions 



There are compelling reasons to think that confinement is a property of QCD, 
and does not result from some other theory. One of such indications is a non-zero 
value of the string tension found in the lattice simulations. As a matter of fact lattice 
simulations give access to many non-perturbative quantities. We are particularly 
interested in knowing the Landau gauge Green functions at low momenta, i.e. at 
energies of order and below than Aqcd- No free quarks or gluons exist at very small 
momenta because of confinement. So, a study of gluonic correlation functions in the 
deep-infrared domain may seem useless. However, in order to study the property 
of confinement from the first principles one has to understand the change in beha- 
viour of Green functions found at low momenta. Knowing these functions exactly 
would be a great support for the future development of the theory, because many 
confinement scenarii (for example the Gribov-Zwanziger scenario) give predictions 
for low energies momentum dependence of the Green functions in Landau gauge. 
Lattice simulations allow to test these predictions. 

Lattice results for different Green functions of QCD have been successfully tested 
at large momenta by perturbation theory up to NNNLO (see chapter|3]and EH/EH). 
We shall see below that lattice Green functions also satisfy the complete ghost Schwi- 
nger-Dyson equation (see Figure l4~2)>. Thus lattice approach gives consistent results 
non only in the ultraviolet domain but also in the infrared one. Of course numeri- 
cal methods could never give us complete Green functions for all possible values of 
momenta. Nevertheless, lattice gives a quasi-unique method for testing different analyti- 
cal approaches, like study of truncated system of Schwinger-Dyson equations, renor- 
malisation group flow equations and other non-perturbative relations like Slavnov- 
Taylor identities. 

Most of analytical predictions are done for the infrared exponents ocp and ocq that 
describe power-law deviations from free propagators when p — > 




(4.1) 



uu 
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where \q p are some fixed parameters of dimension one. When p is large, the func- 
tions F(p ) and G(p 2 ) are logarithmic functions of momentum, and thus a.p g = 0- 
But at low momenta it may not be true. In fact, power-law behaviour is the cru- 
dest approximation, allowing to exhibit the most general features of the momentum 
dependence of the Green functions in the infrared. The real law governing the in- 
frared gluodynamics might be much more complicated. In the following section we 
review (very briefly) different predictions for the exponents ocq and a.p of, respecti- 
vely, the gluon and ghost two-point functions. Next we present our analysis of the 
Slavnov-Taylor identities imposing some limits on these exponents. After this we 
turn to the study of the ghost Schwinger-Dyson equation and test the widely accep- 
ted relation (14.201) between the exponents &q and ocp. Our conclusion (supported by 
numerical simulations) is that this relation is not valid. We revisit the usual proof 
and conclude that either the ghost-gluon vertex behaves unexpectedly in the infra- 
red, or that dp = 0. At the end of the chapter we discuss the results of direct fits of 
two-point functions, and compare our results with other lattice collaborations. 



4.1 Review of today's analytical results 

In this section we quote the main analytical results regarding the infrared expo- 
nents (|4.1|> . We start with the Zwanziger's prediction obtained for gluonic correla- 
tion functions. Next we present the results of other analytical approaches. 



4.1.1 Zwanziger's prediction 



Zwanziger suggested in |50],[5li that the (Landau gauge) gluon propagator va- 
nishes at zero momentum in the infinite volume limit. The argument is the follo- 
wing. The Faddeev-Popov operator (11.681) is positive definite at a local minimum of 
the functional l|1.60|> 

(4.2) 



w,Mp%[U}w) > 0. 



Choosing a test vector 



co a (x) 



.2ne,i 

expi-^x 

A. i 



V 



(4.3) 



where the normalised colour vector % a is an eigenfunction of the "angular momen- 
tum" operator (J b ) a = if ahc , one obtains from dl.681) and l|4.2|) the following limit on 
the mean colour spin 



V 



< 2 tan 



71 

L' 



(4.4) 



Introducing an external colour field source (independent of x), one obtains from 
(H3J) an estimate for the generating functional 



Z(H) 
Z(0) 



z 



— f [DA] e- s[A] +H ?> L * A °> w < e 2V E '' I H '< I tan i . (4.5) 
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The free energy density w(H) = y log Z[H] is convex and bounded from below 
(w(0) = 0), thus one has 

< H7(H) <2£ |H ?( | tan ^. (4.6) 

All connected gluonic Green functions can be obtained by calculating the variations 
of the free energy with respect to the external sources H£(x). The last inequality sug- 
gest that in the infinite volume limit all Green functions vanish at zero momentum. 
However, the inversion of derivation and thermodynamic limit is not supported by 
a rigorous proof. 

4.1.2 Study of truncated SD and ERG equations 

Diverse analytical approaches (study of truncated Schwinger-Dyson equations 
and of renormalisation group equation, see Table l4~T|) agree that the infrared diver- 
gence of the ghost propagator is enhanced, i.e. ocp < ; while they predict different 
values for dlq, mostly around olq ps 1.2. This means that the gluon propagator is 
suppressed in the infrared. However, some groups obtain ocq < 1, i.e. an infrared- 
divergent gluon propagator. Lattice simulations confirmed the prediction for the 
ghost propagator, whereas the lattice gluon propagator seems to remain finite and 
non-zero in the infrared, i.e. ccq = 1. We discuss this question in details in the sec- 
tion 121 



Reference 


Method 






Zwanziger 1 50 1 


see subsectionl4.1.1l 


> 1 


no 


Bloch (13 


SD truncation + perturbation theory 


[0.34,1.06] 


[-0.53, -0.17] 


von Smekal et al. Il53l 


SD truncation 


1.84 


-0.92 


Zwanziger H54I 


SD truncation + Zwanziger condition 


2 or 1.19 


-1 or -0.595 


Aguilar et al. 1 55 1 


SD equation 


0.98 


-0.04 


Kato |56| 


ERGE 


0.292 


-0.146 


Pawlowski et al. 1 57 1 


ERGE 


1.19 


-0.595 


Fischer et al. E5I 


ERGE 


1.02 


-0.52 



TAB. 4.1 — Summary of various analytical predictions 



4.2 Constraints on the infrared exponents and the Slavnov- 
Taylor identity 

Let us consider the Slavnov-Taylor identity ljl.45|> relating the three-gluon vertex 
T\^ v , the ghost-gluon vertex T\u {f> r ) an d the ghost and gluon propagators : 

F(p 2 ) ~ F(p 2 ) ~ 

p A T A}lv (p,q,r) = -^^{5x v r 2 -r x r v )T X}l {r r p-q) - -^^{5 Xv q 2 - q x q u )T Xv {q,p-r). 

(4.7) 



^ 1 1 c i ^ ' i L i ^ rr — 



Taking the limit r — > keeping ^ and p finite, and using the parametrisation G(r 2 ) ~ 
(r 2 ) * G valid for r 2 <C Aq CD , one finds the following limits on the infrared exponents 



dlq < 1 gluon propagator diverges in the infrared, and 

ocp < the divergence of the ghost propagator is unchanged or enhanced in the infrared 

(4.8) 



Let us discuss in details the origin of the limits <4.8D . The ghost-gluon vertex Tu V (p, k; q) 
may be parametrised [59J in the most general way as 



ff c (p,k;q) = f aDC (-iPv)goT V}i (p,k; q ) 



cube 



(4.9) 



/ fl c (-ip v )go- [$v F a(p,k;q) - q v k F b(p,k;q) + p v q fl c(p,k;q) + 



+q v p v d{p,k;q) + p v p v e(p,k;q)] 



(4.10) 



We recall that in this formula — p is the momentum of the outgoing ghost, k is the 
momentum of the incoming one and q = —p — k the momentum of the gluon (all 
momenta are taken as entering). For some particular kinematic configurations we 
use the following dense notations 



a 3(p 2 ) = a(-p,p;0) 

ai{p 2 ) = a(0,-p;p), h(p 2 ) = b(0, -p; p), d x (p 2 ) = d(0, -p;p). 



(4.11) 



The limit r 2 — > leads to an asymmetric kinematic configuration for the three-gluon 
vertex in the l.h.s. of (|4.7|> . This particular configuration allows a general parametri- 
sation Il34l 



T wp (p, -p,0) = (2S pv p p - 5 }ip p v - S pv p }l ) Ti(p 2 ) 



) p p T2{p 2 ) + p F p v p p T 3 (p 2 ). 

(4.12) 

with functions 1^2,3 (p 2 ). The scalar function 7i(p 2 ) is proportional to the gauge 

coupling in the MOM renormalisation scheme. Now, exhibiting the dominant part 
of each term in (|4.7Jl . we obtain 



f((?+r 



Tlil 2 ) (Wv ~ q 2 5iiv) + q 2 T 3 (q 2 )q } iqv + Vip/far) = 
{^(q 2 ) +r 1 (q,r)) (5 }lv r 2 - r }l r v ) + (h{q 2 ) + r 2 (q,r)) q }l (r 2 q v - (q ■ r)r v ) + 



+ (h(q 2 ) + di{q 2 ) + r^qj^r^qy - (q ■ r)r v 



+ 



+ G{q 2 ) 

with ri 2,3 and r\\^ satisfying 



^{(\ 2 ){(\}i(\v - <? 2 V) + tj2 F v(q,r) 



(4.13) 



limri(^,r) = lim^^r) = lim r^(q,r) = 

r^O r^O r^O 

\imi]iuv{q,r) = ]imr]2u V {q,r) =0 

r— >0 r— >0 



(4.14) 



Identifying the leading terms of the scalar factors multiplying the tensors q^q v and 
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{flutfv — q 2 $uv) w e obtain the usual relations ( 1341 ) : 



G(?rw > (4.i5) 

Tsif) =0. 



Using these relations in J4.13|) we get 

Hp 2 ) 



r^o G(r l ) 



a l(q 2 ) (r 2 5 jlv - r H r v ^j + b 1 (q 2 )(r 2 q }l q v - (r ■ q)q }l r v ) = 0. (4.16) 



Thus one sees that if 

ai (q 2 )^0 or b 1 ^0 (4.17) 
then Q4.7JI can only be compatible with the parametrisation l|4.1|l if 

a G < 1. (4.18) 

The condition ()4.17j> is satisfied because at large momentum one has to all orders 
a i(p 2 ) = 1 ( because of the non-renormalisation theorem 181, 11341 ). 

We can also, instead of letting r — > 0, take the limit p — > of ()1.45|l as is done in 
[34J. The dominant part of the l.h.s. of ()1.45|l is 



(28 vv {p ■ q) - p F q v - p v q v ) T^q 2 ) - (s^ - (p ■ q)T 2 (q 2 ) + (p • q)q H q v T 3 (q 2 ) 

(4.19) 

The r.h.s. is the product of F(p 2 ) with an expression of at least first order in p. T\ and 
T2 being different from zero we can conclude in this case that ocp < 0. 

Let us repeat here that all these considerations are valid only if all scalar factors 
of the ghost-ghost-gluon and three-gluons vertices are regular functions when one 
momentum goes to zero while the others remain finite. Under those quite reasonable 
hypotheses one obtains important constraints on the gluon and ghost propagators - 
namely that they are both divergent in the zero momentum limit, and the divergence 
of the ghost propagator is enhanced. 

Let us stress that the limit ()4.8|l on kq disagrees with many other analytical pre- 
dictions quoted in the section l4~Tl 



4.3 Relation between the infrared exponents 

The Schwinger-Dyson equation for the two-point correlation function (and for 
the quark propagator, but we consider only pure Yang-Mills case here) has the sim- 
plest form among other non-perturbative relations between Green functions. It has 
been used to constrain the the infrared exponents. Even more, there is a a commonly 
accepted relation between the infrared exponents 



lap + &g — 0. 



(4.20) 



which we shall discus now. The origin of this relation is the dimensional analysis of 
the Schwinger-Dyson equation for the ghost propagator 



1 



d\ ( F{q 2 )G((q-k) : 



(k-q) 2 -k 2 q 7 
k 2 {q-k) 2 



H^k)), (4.21) 



q 2 (q-k) 2 

where H\(q,k) is one of the scalar functions defining the ghost-gluon vertex : 

q v >T v > v {-q,k;q - k) = q v Hi(q,k) + (q - k) v H 2 {q,k), (4.22) 



where Hi 2 are functions of the factors a,b,c,d,e ()4.9|l . The large momentum beha- 
viour ([34J,[8J) of this vertex depends on the kinematic configuration : 



^■Tf v (-p,0;p) = 1 toaZZorders 



(4.23) 



F-nS(-P/P;0) = i + !fe«s(^) + 



Note that in the case of the vanishing momentum of the out-going ghost (and only 
in this case) the non-renormalisation theorem is applicable |8| and hence 



H 1 (q,0) + H 2 (q,0)=l. 



(4.24) 



Let us now consider two infrared scales k\ = k and k 2 = xk. Calculating the diffe- 
rence of the Schwinger-Dyson equation ()4.21|> taken at scales k\ and k 2 and suppo- 
sing for the moment that a.p 7^ one obtains 



F(k) F(Kk) 



a (l-K- 2a ^)(k 2 )-^ 



d\ (F(q 2 ) f(k-q) 2 2 



(In) 



k 2 



x 



G{{q - fc) 2 )Hifa,*) G((q - Yk) 2 )H^q,Kk) 



((q-Kk) 2 ) 



2^2 



(4.25) 



This integral equation, as well as the initial equation l)4.21|l , is written in terms of 
bare Green functions, and the integral may contain ultraviolet divergences. It can be 
cast into a well-defined renormalised form by multiplying (in ()1.38|l ) (resp. F^) 
by Z^ 1 (resp. Z^ 1 ) and the bare coupling by Z g 2 = Z3Z 2 , and finally multiplying 

the k 2 term by Z3. However, in the subtracted equation ()4.25|l all ultraviolet diver- 
gences are cancelled, as well as the Z3A: 2 term. Thus the subtracted Schwinger-Dyson 
equation holds both in terms of bare and renormalised Green functions without any 
explicit renormalisation factors. 

We now make the hypothesis that there exists a scale qo below which the power- 
law parametrisation is valid 



G(q 7 



2\a G 



Ftf 



[q 2 r, for q 2 <q 2 . 



(4.26) 



The equation l)4.24|> suggests that if both functions Hi 2 are non-singular then one 
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can suppose Hi(q,k) ~ 1 in H4.25I , and H4.2QI> is straightforward by a dimensional 
analysis. However, we have a priori no reason to think that the scalar functions 
Hi (q, k) and Hiifl, k) are separately non-singular for all q, k. Writing for example 1 

Hi fa,*) ~ if)*' h ( q -£,K) , (4.27) 

with a non-singular function h\, we keep all the generality of the argument admit- 
ting a singular behaviour of the scalar factor H\(q,k). Doing the dimensional ana- 
lysis of the equation (14.25b without putting H\(q,k) ~ 1, we obtain that the relation 
(14.20b is satisfied if and only if the following triple condition is verified II60II : 

2cc F + a G = < — ► < dp = (4.28) 

[ dp + kg < 1 

All possible cases and limits obtained from the integral convergence conditions are 
given in Table l4~2l As we shall see the case 2 is excluded by lattice simulations. The 
case 4 is particularly interesting, it corresponds to the situation when the power- 
law infrared behaviour of the ghost propagator is the same as in the free case, and 
no relation between the infrared exponents follows from the Schwinger - Dyson 
equation. We shall return to this the discussion of this case in the section l4~5l 



case 


1 

0L F £ 

«F + occ + «r < 1 


2 

a F / 

OCp + OLq + Oil > 1 


3 

dp = 
#G + a r < 1 


4 

«f = 
OiQ + OiY > 1 


IR 










convergence 


Up + &y > —2 


Dip + Oil > —2 


ar > -2 


dp > -2 


conditions 


«G + a r > —1 


Uq + Oil > —1 


ftic + #r > —1 


etc + a r > —1 


SD 

constraints 


2&F + OCq + DCy = 


dp = —1 


excluded 


none 



TAB. 4.2 — Summary of the various cases regarding the a coefficients 



The first and the last conditions (14.28b are compatible with limits coming from 
the analysis of the Slavnon-Taylor identity ()4.8|) , and are also consistent with lattice 
simulations (see section 1431 Il60l0 . If one of the conditions (|4.28[> is not verified then, 
according to the Table 14.21 (14.201 1 should be replaced by 

lap + Di G + dp = 0. (4.29) 

In the following section we present the results of a numerical test of the relation 
(14.20b , and thus we probe the validity of the condition on ctp. 

One remark regarding the power-law parametrisation is in order. Suppose for 
the moment that this parametrisation is exact below the scale qo defined in (14.26b . 
Then one can differentiate (14.25b n times with respect to k, keeping q, k finite. We 



fact there are many possible parametrisation. We choose (14.271 in order to illustrate the argu- 
ment that follows. 



WW 



V_ 1 Id j^y ill t. ± — 
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obtain 

*) ^ (_2«p) ■ . . . ■ (-2a F - n) cc / ^ f ^-k*) 2 )^) 

(4.30) 

The r.h.s of the last equation is not equal to zero for finite A:, and thus one immedia- 
tely has 

a F ^-|, n = l,2,-... (4-31) 

Thus any half-integer predictions for ap should be considered as an indication of 
incompleteness of the power-law parametrisation ()4.1|l . 

4.4 Lattice study of the ghost Schwinger-Dyson equa- 
tion 

4.4.1 Complete ghost Schwinger-Dyson equation in the lattice for- 
mulation 

In order to derive the discretized version of the ghost Schwinger-Dyson equation 
we repeat the same steps as in the continuum case dl.321 - 11.361) but for the lattice 
version of the Faddeev-Popov operator J1.68I I. We define the covariant Laplacian 

Ag = (Gj*(x) (S x , y - S yrX+e!l ) - Gf(x - e H ) {s y , x - ep - 8 y , x ) ) . (4.32) 

The appearance of Ajj in ()1.68|l is due to the appropriate discretisation of the usual 
Laplacian operator A, dictated by the non-locality of derivatives in the lattice for- 
mulation, i.e. replacement of the V operator by its covariant version. 
Multiplying (11.681 1 by (x, y) from the right, one obtains 



2(N? - r 



j^*n&*)qs(u;z,x) = s y , u - 



This is an exact mathematical identity for each gauge configuration U, and thus 
the consequences that can be derived from this relation are free of any ambiguity 
originating from the presence of Gribov copies. Performing an averaging (•) over 
the configurations U one gets 

1 - /. , _(2) 



' rabc ( 4 - 34 ) 

- 2(4^1) (^^ F S5(^y + - 4(y - e x) F Z(u>y - 

This averaging procedure depends on the way chosen to treat the Gribov problem : 
the particular set of configurations over which it is performed depends on the pres- 



cription which is adopted (fc/bc procedures on the lattice, restriction to the fun- 
damental modular region ; see the subsection 12.4.21 for details). Consequently, the 
Green functions may vary but they must in any case satisfy the above equation, 
even when the volume of the lattice is finite. 

Like in the continuum case, we perform a Fourier transform and obtain : 

^TrXV>* (A u (0 /Z )F 1 2 nf (U, Z ,x)) = 1 - ifafo)-J^ (^(OJfg&lI,?)) 

(4.35) 

Although the equations (14.331) and (14.341 1 have to be exactly verified by lattice data, 
the relation (J4.35|> does only approximately (within statistical errors) since it relies 
on translational invariance, which could be guaranteed only if we used an infinite 
number of Monte-Carlo configurations. 

The presence of Au in the last equation is due to non-zero lattice spacing effects. 
Indeed, lattice perturbation theory possesses an infinite number of ghost-gluon ver- 
tices depending on the lattice spacing a, giving tadpole contributions like the one 
presented at the Figure l4~Tl Such tadpole contributions maybe estimated by a mean 



2 

■a 



FlG. 4.1 — Example of the terms in the Schwinger-Dyson equation on the lattice. 



field method II6T1 . Using the average plaquette (P) (for f> = 6.0 (P) ~ 0.5937) one 
predicts a tadpole correction factor <x (p)~(i/ 4 ) ~ 1.14. These terms disappear in 
the continuum limit, but they do so only very slowly : the tadpole corrections (1 - 
plaquette) vanish only as an inverse logarithm with the lattice spacing. This is to be 
contrasted with the corrections arising in the r.h.s which are expected to be of order 
a 2 . Our lattice calculation [60 J gives 

A u ~ A/ (1.16 ±0.01), (4.36) 

almost independently of the momentum. This is in good agreement with the correc- 
tion factor 1.14 quoted above. 

We see from Figure 14.21 that the lattice Green functions match pretty well the 
SD equation ()4.35|> in both the ultraviolet and infrared regions. Lattice propagators 
were successfully checked by the perturbation theory at large momentum, and they 
satisfy the ghost Schwinger-Dyson equation. This means that lattice approach gives 
consistent results also in the infrared. 
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FlG. 4.2 — Checking that lattice Green functions satisfy the ghost SD equation \4.35\ . We plot jj^F (p 2 ) 
gojgtlf^iA^O) Fg$(A, p)) compared to 1. 



4.4.2 Checking the validity of the tree-level approximation for the 
ghost-gluon vertex 

The simplest approximation (used by many authors, see section l4~T|> of the ghost 
Schwinger - Dyson equation l)4.21jl corresponds to the case 

H 1 (q f k)=l Wq,k, (4.37) 



an approximation motivated by the non-renormalisation theorem l)4.24j> valid for the 
sum Hi (q, k) + Hj{c\, k) when k = 0. This gives 

J_ -i | jNc f dS f F(q 2 )G((k-q) 2 )(k.q) 2 -k 2 q 2 x 

F(k) k 2 J (2tt) 4 V q 2 {k-q) 2 (q - k) 2 )' K ' 

Strictly speaking this equation, written in this way, is meaningless since it involves 
UV-divergent quantities. However it is well defined at fixed ultraviolet cut-off. 

We want to check whether lattice propagators satisfy it. According to pertur- 
bation theory, it should be approximately true at large k. Lattice propagators are 
discrete functions of momentum and thus one has to handle the problem of the nu- 
merical evaluation of the loop integral I in (14.381 1 . Let us express the integrand solely 
in terms of q 2 and (k — q) 2 



I 



d 4 q F 2 (q)G(k-q) l ( k-q) 2 | (k 2 ) 2 + (q 2 ) 2 - 2k 2 q 2 q 2 + k 2 



(2tt)4 q\k-q) 2 
Then we write 



4(k-q) 2 

I = h + I 2 + 1 3 + h + k + h, 



(4.39) 



(4.40) 
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each corresponds to one term in (14.39ft . All these integrals have the form 

Ii = Ci(k) I A^f^hiik-q). (4.41) 

The convolution in the r.h.s. is just the Fourier transform of the product at the same 
point in configuration space : 

/ -0ff^^ k ~^ =F + (F-(/i-)MF-(^)[x])(fc), (4.42) 

where F_ (/) (x) is an inverse and F+(/) (k) a direct Fourier transform. Thus, in order 
to calculate the integral I from discrete lattice propagators one proceeds as follows : 

1. calculate {/;}(p) and {hi} (p) as functions of F(p), G(p), p 2 for all i 

2. apply the inverse Fourier transform F to all these functions and get fi(x) and 

h{x) 

3. compute the product at the same point fi(x) ■ hi(x) 

4. apply the direct Fourier transform F + to fi(x) ■ hj(x) 

The integrands in |J4.39|) depend only on the squared norms q 2 and (k — q) 2 , and 
thus the angular part may be integrated out, giving the four-dimensional Hankel 
transformation 

i r°° 

/(W)[P]=A / h(\p\r)r 2 f(r)dr 
\V\ Jo 

/W = (2^7 i h{pr)p 2 f{\x\)[p]dp. (4.43) 
These integrals are evaluated numerically by means of the Riemann sum 

fir) = (2tt) -'M- 1 £ h(r Pi )p ^ [Pi] (Pi ~ Pi-i)' ^o = 0, (4.44) 

i=i L 

where N is the number of data points. The inverse transformation is done in the 
similar way. In practice, because of the lattice artifacts (see subsection I2.3.2|l which 
become important at large p the summation has to be restricted to p < p max — 2.2 
instead of the maximal value 2n. 

Now we are ready to check the approximate equation (14 .3811 on the lattice. We 
still have to face the same problem we have already encountered in the previous 
subsection, namely that the lattice Faddeev-Popov operator involves the non trivial 
discretisation of the Laplacian operator. This is taken into account by means 
of the substitution of A u (p 2 )/p 2 to the "1" term in the l.h.s of equation d4.38|) We 
present on Figure 14.31 the result of the numerical integration described above. We 
have chosen for this purpose the data set from the simulation with the gauge group 
SU(3) at £ = 6.4, V = 32 4 ,a~ 1 w 3.6 GeV. One sees that the equality is achieved at 
large momenta, but in the infrared the naive approximation of the ghost Schwinger- 
Dyson equation fails. The errors on Figure l4~3l include statistical Monte-Carlo errors 
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FlG. 4.3 — Checking whether lattice Green functions satisfy the ghost SD equation \4.21\ with an assumption 
H\(q,k) = 1. The upper line(circles) correspond to the loop integral in \4.21\ . and the down line (triangles) 
corresponds to 1/F(p 2 ) — 1. In this plot a -1 « 3.6 GeV. 

for F(q 2 ) and G(q 2 ) and the bias coming from the UV cut-off of the integral I. 

We see that at small momenta (below w 3 GeV) the ghost Schwinger-Dyson 
equation with the assumption H\(q,k) = 1 is not satisfied. However, it is quite 
difficult to establish whether this disagreement is due to the infrared or ultravio- 
let dependencies of H\(q,k). To check this one has to know H\(q,k) for all values of 
q, k. Unfortunately this information is not available. Thus the main conclusion of the 
present subsection is that the scalar function H\(q, k) plays an important role in the 
infrared gluodynamics, and it cannot be set to one. 
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4.5 Direct fits of infrared exponents 

We have seen in the previous section that lattice simulations give consistent re- 
sults for the Green functions at all momenta. Another interesting feature that we 
have established is the important role of the scalar factor H\{q,k) coming from the 
complete ghost-gluon vertex ()4.22[) . In this section we present numerical results allo- 
wing to check the relation (I4.20I ). After this we present our results for direct fits of the 
exponents ocp and ocq, and compare them to the results of other lattice collaborations. 



SU3 R=5.75 Vol=32 



I 

0,5 




FlG. 4.4 — Direct test of the relation 2ap + ocq = 0. If the last is true F 2 G has to be constant in the infrared. 
We see that it is clearly not the case. In these plots a -1 k 1.2 GeV, so the peak is located at m 600 MeV. 



4.5.1 Testing the relation 2ctp + olq = 0. 

In order to test the relation (14.201) we plot at Figure l4~4l the quantity F 2 (p 2 )G(p 2 ). 
If all the conditions ()4.28|> are satisfied this quantity should be constant in the infra- 
red (or slightly varying). We see from Figure l4~4l that in the infrared (below rs 600 
MeV) the quantity F 2 G is not constant, and thus one of the conditions (14.281) is not 
verified. We have seen that the conditions ocp ^ and ocp + Kq < 1 are consistent 
with the limits (14.81 ) from the Slavnov-Taylor identity l|4.7|). We have also seen (cf. 
Figures l4~2l and l4~3)) that neglecting the momentum dependence of the vertex is not 
possible in the infrared, because in this case the ghost Schwinger-Dyson equation is 
no longer satisfied by lattice propagators. Thus the only possibility is to admit that 
Hi (q, k) plays an important role, and that the relation (14.201 ) is not verified. If dp ^ 
then the modified form (14.291 ) that takes in account the singularity of Hi (q, k) should 
be considered (according to Fig J4.3|) , with dp < in our parametrisation. This sin- 
gularity is probably related to the non-perturbative power corrections to the vertex 
discussed in the subsection 13.2.31 

Another reason to think that the relation (14.201 ) is not exact is the dependence of 
ocp and o.q on the choice of the Gribov copy. We have seen in the section l2~4l that the 
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• SU2 [3=2.3 Vol=48 4 

■ SU2 P=2.3 Vol=32 4 

Vol=48 4 Fit G <2) =(p 2 ) ao "'(a+bp 2 ) a Q =1.004±0.015 

— - Vol=32 4 Fit G <2) =(p 2 ) v '(a+bp 2 ) a G =0.968+0.01 
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FIG. 4.5 - G^(p 2 )from lattice simulation for SU{2) (left). f> SU (2) = 2.3 and p S u(3) = 5.75. The volumes 
are 32 4 and 48 4 for SU(2). In these plots a' 1 w 1.2 GeV. 



low-momentum dependence of the gluon propagator is not sensitive to the bc/fc 
choice while the infrared behaviour of the ghost propagator depends on it. But the 
Schwinger - Dyson equation for the ghost propagator is independent of the choice 
of the copy, because it is valid exactly for every gauge configuration, even on a fi- 
nite lattice (see equations dl.361) and (I4.34D V Hence if there is a relation between the 
infrared exponents ocp and ocq resulting from the ghost Schwinger - Dyson equation 
then it could not depend on the choice of the copy. Thus it is not possible to have a 
relation with ocp and o.q alone. This above argument is not directly applicable in the 
case ap = 0. 

According to the analysis performed in the section l4~3l and given (see next sub- 
section) that the case 2 of the Table 14.21 is excluded by lattice simulations the follo- 
wing explanations of the non-validity of the relation 2&p + &q = are possible : 

1. The ghost-gluon vertex contains scalar factors that are singular in the infrared, 
i.e. dp 7^ in the equation 14.271 

2. The case 4 of the Table 14.21 is realised 11621 and hence there exists no relation 
between the infrared exponents. Let us recall that in the above case one has 
tx,p = and k.q + dp > 1. If the ghost-gluon vertex is regular in the infrared 
then one has 




(4.45) 



lAtO \J1 11 111 CU. C CA|/U1 LCA LLO 



4.5.2 Lattice fits for otp and dlq. 

Let us now discuss the direct fits for the infrared exponents ocq and dp. The 
examples of such fits of lattice data are presented on Figure l4~5l The errors are quite 
large, leading to an instability in the fit results. That is why we fit both propagators 
in the infrared to the formula 

( q 2 n\ + M 2 ) (4.46) 

where we added an additional term of the form jiq 1 in order to describe a situation 
like the one at Figure l4~3l lef t) where G^ 2 \p 2 ) seems to go to a finite limit when p 
goes to zero. The obtained values for ocpQ are summarised in Table l4~3l For SU(2) 



Group 


Volume 






OCp 


SU(2) 


48 4 


2.3 


1.004(15) 


-0.087(15) 


SU(2) 


32 4 


2.3 


0.968(11) 


-0.109(14) 



TAB. 4.3 - Summary of the fit results for the F and G functions 

and the larger lattice volume the value obtained for ocq is compatible with l.We 
also take into account our experience from previous studies of the gluon propagator 
where we have always observed that the gluon propagator goes continuously to a 
finite limit in the infrared region (see Figure l4~6b . However, the fits are quite instable, 
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FlG. 4.6 — The continuity of the lattice gluon propagator in the infrared. 

and depend a lot on the choice of the fit formula that can considerably change the 
result. The main problem is the lack of data points at low momenta. 

Regarding the gluon propagator another strategy may be taken. It consist in ex- 
trapolating the available data to the infinite volume limit. A very detailed study of 
the gluon dressing function and specially of its volume dependence at k = has 



/ -r 



already been performed in 1631 . This study shows that a value ocq = 1 is compatible 
with the data (the dressing function shows no signal of discontinuity in the neigh- 
bourhood of zero) and that no pathology shows up as the volume goes to infinity 
Let us compare all available lattice results for the point G^(0) and check whether 
there is an agreement between the data. Following [63 1 we renormalise the gluon 
propagator in the MOM scheme at 4 GeV and use the suggested extrapolation for- 
mula 

Gf (0,p = 4GeV) = G^l(0 /F = 4GeV) + (4.47) 

We compare the results of l63l . l64l and our data from the Table l4~4l The results for 
the fit parameters (0) and c are presented in the Table l4.5l We are aware that not 





V in units of a 


bare propagator G^(p) 


1/aL in GeV 


5.7 


16 4 


16.81 ±0.13 


0.0672 


5.7 


24 4 


15.06 ± 0.29 


0.0448 


5.8 


16 4 


19.12 ±0.16 


0.0841 


5.9 


24 4 


18.12 ±0.30 


0.0685 


6.0 


32 4 


17.70 ± 0.59 


0.0615 


6.0 


24 4 


19.67 ±0.35 


0.0821 



TAB. 4.4 — Physical lattice sizes and raw data for the gluon propagator at zero momentum (p) from our 
old data. 



reference 


G^ 2) (0, y. = 4 GeV) in GeV 2 


c in GeV 2 fm 4 


max vol in fm 4 


L63J 


7.95 ± 0.13 


245 ± 22 


2000 


Table El 


9.1 ±0.3 


140 ± 50 


90 


L64J 


10.9-11.3 


47-65 


110s 



TAB. 4.5 — Summary of the infinite volume zero momentum propagator and its slope in terms of \ IV for 
three different simulations. The largest volume used in the fit is also indicated. The statistical error is not quoted 
in H. 



all systematic errors are taken into account : O(a) effects, effect due to different lat- 
tice shape, insufficiently large volumes (for the second and third lines), uncertainty 
in the estimate of the lattice spacing in physical units, etc. However, it seems that 
not only there is a clear indication in favour of a finite non vanishing zero momen- 
tum gluon propagator, but that different lattice collaborations agree on the value. Of 
course a more extensive study is necessary to check this statement. The other free 
parameter of the fit - the slope c - is clearly different, but still all the values are in 
agreement in the order of magnitude. 

We conclude [651 thus that all available numerical results point towards a finite 
non-vanishing and zero momentum renormalised lattice gluon propagator in the 
infinite volume limit. This suggests that ocq = 1. An additional study at much larger 
lattices is needed to get a reliable result for this infrared exponent. 

This last result is in conflict with the limits found from the study of the Slavnov 
- Taylor identity (14.81 1, and contradicts the Zwanzigers's prediction that the gluon 



propagator is infrared suppressed. However, it is very close to the results presented 
in Il55i 

Let us finally discuss the gauge-dependence of the parameter ocq. This is still an 
open question. However, the results of works |66J,|67J suggest that the value of ocq 
does not change drastically (see Figure 14.8b when changing the gauge parameter £. 
This question deserves a separate study. 




p, GeV 



FlG. 4.7 — The fit of the ghost scalar factor F(p) to the formula A + Blogp. It suggests that the infrared 
divergence ofF(p) is very slow. Hence ocp is close to zero \62\ . The simulations was performed on aV = 32 4 
lattice at f> = 5.8. The Landau gauge was fixed using the fx. choice for the Gribov copies. 

To finish this chapter, let us summarise the lattice results. We have found that 
txp is very close to (see Figure W7\. &q is close to 1 and the widely used relation 
2ocp + ot-G = is not true. Going back to the possibilities given in the Table l4~2l we 
find that the cases 2 and 3 are not realised. We are left with the cases 1 and 4, and 
for the moment lattice simulations cannot say which possibility is true. However, 
all numerical results are better explained by the possibility (14.451 1 corresponding to 
the case 4 of the Table l4~2l supplied with an hypothesis of the regularity of the scalar 
factors entering the ghost-gluon vertex. We recall that in this case one has : 

a F = 

a G > 1. (4.48) 
no relation between ocp and ocq follows from the ghost SD equation. 

Note that it is still in conflict with the constraints coming from the Slavnov - Taylor 
identity H4.8|>. Thus the essential question today is to understand whether ccp = or 
not l62l . And, of course, a study on larger lattices is necessary to perform better fits 
of the infrared exponents. 
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FlG. 4.8 — Transverse part of the gluon propagator p 2 G^ 2 \p) in covariant gauges as a function of p. The 
two sets of data refer to (£ = A) £ = (Landau gauge) and £ = 8, 221 thermalized SU(3) configurations at 
f$ = 6.0 with a volume V x T = 16 3 x 32. Extracted from l66t . 



Chapitre 5 
Conclusions 



In this chapter we discuss the conclusions of the present dissertation. Lattice si- 
mulations is a great tool to study the non-perturbative effects in QCD. The main goal 
of this dissertation is to exhibit how these effects influence the momentum beha- 
viour of different Green functions in Landau gauge. Our hope is that the knowledge 
of the change in momentum behaviour at low momenta can help in the understan- 
ding of one of the most difficult puzzles of QCD - the mechanism of confinement. 

In the chapter[3]the large momentum behaviour of the ghost and the gluon propaga- 
tor of a pure Yang-Mills theory in Landau gauge is investigated. The main parame- 
ter under study is the scale Aqcd- We show that the values of Aqcd fitted from the 
ghost and the gluon propagator are consistent. However, the available momentum 
range (from « 2 to ~ 6.5 GeV) is situated in the zone where non-perturbative effects 
cannot be neglected. So at first glance the agreement between Aqcd extracted from 
different Green function may seem strange. An explanation of this fact comes from 
the OPE analysis allowing to estimate the influence of dominant non-perturbative 
power corrections. We found that these corrections are the same in the case of the 
ghost and gluon propagators, that is why the values of Aqcd are compatible. Accor- 
ding to the OPE calculation, the value of Aqcd that is extracted from the propagators 
is modified by the non-perturbative effects. We used the fact that the equivalence of 
the leading power corrections implies that their ratio is free of power corrections at 
the considered order. Thus the ratio of the ghost and gluon propagators is a quantity 
which is better described by perturbation theory in the considered energy interval 
than the propagators themselves. Indeed, our analysis of the lattice data showed 
that the fit of the ratio gives a smaller value for Aqcd (~ 270 MeV in the MS re- 
normalisation scheme) compared to the value obtained from the separate fits of the 
propagators (A^ ps 330 MeV). Both fits are performed on the same data samples. 
This speaks in favour of a presence of non-perturbative power corrections in the 
interval [2 GeV, 6.5 GeV], in agreement with the OPE predictions and the results of 
previous investigations. 

One can use the Slavnov - Taylor identities in order to find other relations bet- 
ween the Wilson coefficients in different Green functions. For example we found 
that the dominant oc \ power corrections are the same in the case of the three gluon 
vertex and the ghost-gluon vertex in asymmetric kinematic configurations (with the 



gluon momentum set to zero). We have also shown that the power correction to the 
ghost-gluon vertex with vanishing momentum of the entering ghost is equal to zero 
(in Landau gauge). However, this is not true if the momentum of the entering ghost 
is not exactly equal to zero. In this case the vertex has an <x \ power correction that 

becomes quite important at low momenta. 

As a partial conclusion we stress the attention of the reader on the fact that lattice 
simulations are very successfull in describing the Green functions at large momen- 
tum, the results are consistent with the predictions of perturbation theory, comple- 
ted with the OPE calculation of the power corrections, up to NNNLO. 

In the chapter 21 we turn to the study of very low momentum behaviour (of order 
and below Aqcd) of the Green functions. One of the very interesting puzzles in the 
infrared is the problem of the Gribov ambiguity. The lattice method has an advan- 
tage to explicitly perform the Gribov quantisation. However, the gauge is not fixed 
in a unique way and there are Gribov copies on the lattice. We showed that the 
probability to find a secondary Gribov copy possess the property of scaling. This 
probability increases significantly when the physical volume of the lattice excesses 
some critical volume (of around (2.75/ v 7 ^) 4 i n the case of SU(2) gauge theory). Our 
conclusion is that in order to study the non-perturbative effects one has to work at 
lattices with physical volume larger then the critical one that we found. 

Our first step in the study of low-momentum behaviour of Green functions is 
to check that lattice simulations can give reliable results in the infrared. For this 
purpose we verified (numerically, see Figure l4~2)> that lattice Green functions satisfy 
the complete ghost Schwinger - Dyson equation ()1.38|l for all considered momenta. 
These tests allow us to conclude that numerical simulation on the lattice give rele- 
vant results not only in the ultraviolet domain, but also in the infrared one. 

The quantitative parameters we are interested in are the infrared exponents ap 
and ocq describing the power-law deviation (this is a crudest approximation) from 
free propagators (ghost and gluon respectively) in the deep infrared. Our analysis of 
the Slavnov-Taylor identity relating the three-gluon vertex, the ghost-gluon vertex 
and the propagators showed that the power-law infrared divergence of the ghost 
propagator is unchanged or enhanced in the infrared (compared to the free case), 
and that the gluon propagator must diverge in the infrared. The latter limit is in 
conflict with today's lattice results yelding a finite non-zero gluon propagator at 
zero momentum, and with most present analytical estimations that we quote in the 
section |4~T1 that support a vanishing gluon propagator in the infrared. 

Another analytical relation imposing constraints on the infrared exponents is 
the ghost Schwinger-Dyson equation. We revisited the commonly accepted relation 
()4.20J> between these exponents saying that 2a p + ag = 0. According to our ana- 
lysis this relation is true only if the ghost-gluon vertex contains no singularity, in 
none of the scalar functions defining the vertex. Our numerical studies showed that 
the relation in question is not valid, because F 2 G is infrared suppressed, and hence 
2ap + ag > 0. This statement is supported by the fact that lattice propagators do not 
match the reduced Schwinger - Dyson equation (see Figure I4.3|l , whereas the com- 
plete one is perfectly verified, Figure 14.21 There are even more reasons to think that 
the relation (14.20b is not true. First, we have seen that the non-perturbative ghost- 



gluon vertex contains singular contributions from the (A ) condensate for most ki- 
nematic configurations. However, it is difficult to estimate the role of such contribu- 
tions at very low momenta. Second, we have seen that the infrared behaviour of the 
ghost and gluon propagator seem to vary differently with the choice of the Gribov 
copy (be or fc). We have seen that the form of the complete ghost Schwinger - Dyson 
equation does not depend on the choice of the Gribov copy. If ap ^ and it depends 
on the choice of the copy, it is impossible to have an exact relation between ap and 
OiQ alone, and to satisfy the condition of independence of the choice of the copy. 

The direct fit of the propagators in the infrared supports the prediction that the 
infrared behaviour of the ghost propagator is enhanced in the infrared. But this 
enhancement is very slight. The gluon propagator is found to be infrared finite. This 
last result is in conflict with the limit found from the analysis of the Slavnov - Taylor 
identity. For the moment we have no explanation regarding this disagreement. 

Summarising the numerical results, we found that the gluon propagator is finite 
in the infrared {ctQ £s 1), that the infrared divergence of the ghost propagator is 
almost the same as in the free case (ap « 0) and that the commonly accepted relation 
2a p + ag = is not true. Going back to the analysis of the ghost Schwinger - Dyson 
equation (see Table l4l2)l. two solutions are possible : 

1. The infrared exponent of the ghost propagator ap is strictly equal to zero, i.e. 
the power-law infrared dependence is the same as at large momenta. This im- 
plies that no relation between ap and ac follows from the ghost Schwinger - 
Dyson equation. If we now suppose that the ghost-gluon vertex contains no 
(infrared) singular components then all our lattice results are perfectly descri- 
bed. 

2. The infrared exponent of the ghost propagator ap ^ 0, then there is a relation 
between ap, aQ and ap following from the ghost Schwinger - Dyson equation. 
The fact that the relation 2ap + ac = is not verified on the lattice suggests 
that there is a singularity in one of the scalar factors defining the ghost-gluon 
vertex i.e. ap ^ 0. 

The today's lattice results speak in favour of the first possibility, but calculations at 
larger lattices are necessary in order to conclude. 
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